





LSI: 35) Yayo -F 


RANGELAND 
MONITORING 


x } WAC 












.t A; 









































Statistical Considerations 


@ py 
yi 


TR 4400-8 ¢ 1992 


J 














United States Department of the Interior 
Bureau of Land Management 


Statistical Considerations 
in 
Rangeland Monitoring 


by 


Frederick K. Martinsen 
BLM Service Center 


and 


John Willoughby 
California State Office 


Technical Reference 4400-8 
September 1992 


BLM/SC/PT-92/010+4400 











Table of Contents 

















Section Page 
- i TEETER Snr ne hie e ates inst tee nena NI l 
TT NS TE TET a Te Ie TT TD Te TNO 3 
ee 3 
EL Te eT a a 3 
i ala eraclnrentnal 4 
TEETER ee eT en TT ET en IEE: 5 
ee nt ITI TINIE... ccccncnncesensstennssonensnsnnsnenenenesnssetsssoonessenssossneee 7 
EE 8 
rw FT |_| TEE 9 
2.8 Where and how should the information collected be stored? ................:.sccsseseeeeeeeees 10 
Tg Tee EE TE 11 
3.1 Methods for photo plot, line intercept, and step point ...................scsseeseeseseeeeeeseeeeeeees il 
| TTS Se a ee Tne EEE 11 
Ge = ETT TE ne EEN 13 
 - Si EERE ES arene mea L ET eT ne eT TERETE EEE EE 13 
LETTE ane ae nT REET ETE 13 
| | a — | nnn Ne eee ETE E 13 
3.1.6 Analysis of one transect or plot per year ................cccccccceesecesseeeseseeeeseeseeeeeeeeenees 14 
3.1.7 Analysis of two or more transects Per Yea ...............ccccceessesessesseeseeeeeseeseeeeeeeeeees 14 
3.2 Methods for community structure analysis and Daubenmiire ......................::00seeeeeeee- 14 
ee TIT sccncsnscccnntenstnnindmnenenatienteenesnesennrtnhenessnesestanassecnenaneecennnansenes 15 
ee II ccncscnterecccnssesns -cxnsnssennsssesnensnsensenincesseescscnecencescenssecasocnteene 16 
Ee OIE a cecccnscncsscenenscinsnecevccnssescensnneantnnesecssneoesesentensnesonsnsenensees 16 
Be | OI ccissnscasnscncnncnnsessenseunsnsnsesnsesosstnnenetanecescneensessessemtesonenessenes 16 
Tee GiURIIIIEIIIIIIIIIIIE ¢00/cnncnsstnssusenenindntensnsenssneseenansohnepeoonsenmnecnnsessorennsentenstonane 17 
ee I Ge Fe GO ccccccesevcccccceseccesscvscecesecsesocsccsccoconscescncnscscosenscnesoecoesioees 17 
III os ncrtiensenipnonneconsnnntnenscesennensnsnenseanneuneenerensseennscessoonemmnennnntes 17 
3.2.8 Analysis of one transect Per Year ............cccsccssesseseeseeeeeseesessssessesssseesesesssseeeseeeesens 18 
3.2.9 Analysis of two or more transects Per YOAF ............cccccceseeeseeseseeseseesesseseensceneeneeees 18 


+e 














3.3 Methods for pace frequency, quadrat frequency, and nested frequency ................. 19 

















5 anne aC aan N EEE 19 
I a cacnalall 19 
I 20 
I i ll 20 
| TTT ITT a ae TE ETE TEER 21 
3.3.6 Analysis of two or more transects per Year .................scesceescesseseeeseeeeeseeseeseeseeneeneees 21 

©, Ce Ce CI I oecccccccccccececcncccccnesecceccescrsccceveccsncocececcccenesseceecsooesees 22 
TT Ea 22 
4.2 Data preparation in all other methods ......................:cscsscsessssecseceeeeesecseseeeseseeceneeeees 23 
LTT ETAT TET aT TENET TEE 23 

— ELST AE en 23 
ER EA on Oe a ET 23 

SF | EET nen ene een SE 24 
at elmnnnnteinetnnmniannentenntinitel 24 
FT na ane ee Ee .24 

5.3 Example of Two-way ANOVA (Randomized block design) .................:.csseeeeees 25 

5.4 Example of Duncan’s Multiple-Range Test .....................csccsccsseseseseeesereeeeeeeeeeeeeeees 26 

5.5 Example of Two-way ANOVA with repeated measurements ......................sc000000+ 29 

5.6 Example of Two-way ANOVA with more than one observation per block ........... 29 

ee a eiiiisiinestiscatenscnsciisiinnsenenensnstnieidaiiiinininniaenamiatannninnnniiimimianmsimanions 30 

SE earanTnTTUIT UUITTTTTTTTT Tis iin ciineainr <eeliheheeiadtntalieidanenmnnniniitadiintatias 31 
Er TTIT TIED aiitticsiethaiaiacipantienecrernimencesicatbiniuasidaieaeiensiciguaiaineiennieeeameinnniinitendanetaitemmmnnneceaeennantiemnieniinte 33 
SE deininetennnietprtnnneensisienanestanennenesnnan cennemninnisiiamentaentnanemenninnnnennninnamnnnnn 35 
Re ICID ccnicristieicnsnbncnnrnmanninatencemesnenneennnesinanintennneiinesnmmmameetint E-36 
Pace ies ticrennessctnesbesinanninmnnitiaperenimnennenteenennemsenieninntieimnntmemmnensnenneninanantints E-36 
IIE einstenncansusanensnsicnnstsnnnsrenntnereneennseneesnuseneeseneensennstenssiennenanensnenennnnntani E-37 
E,W RRIIIIIIIIGIED cocoscsssnusssensnossneencnssscevensessttnnsnnsneunnnsoccseessssenssesestentenenenteqeocnscnsestentiets E-38 
Sh, FITTING ccnsnsscisinsnistitenenenngeinssnemusinnetimnatenenmeesnenenetnnnssennessennetarsanantmiancenacetettt E-38 
6. SPSSX command file for t-test with independent samples ...................:..0ccseseeeeeeeseeees E-39 
7. SPSSX results for t-test with independent samples ....................:.:c0sccssseesseeereeeneeeeees E-40 
8. SPSSX command file for t-test with paired samples ...............cccccssseseereeerseeseeeeeenenees E-41 
9. SPSSX results for t-test with paired samples .................:csseseeeereeeeeereseeeeeeeseneesees E-42 
10. STPK results for t-tests with independent and paired samples...................c.sceseeseeeees E-43 
11. SPSSX command file for Two-way ANOVA with randomized block design............ E-47 
12. SPSSX results for Two-way ANOVA with randomized block design .................0006 E-48 

i 











13. STPK results for Two-way ANOVA with randomized block design ......................... E-50 
14. SPSSX command file for Duncan’s Multiple-Range Test .....................----:c-sesseeeeeeeees E-54 
15. SPSSX results for Duncan’s Multiple-Range Test ........................-ccccseeeeseeeeeeeeeeeeeeeeees E-55 
16. /RENO/XMONITOR results for Duncan’s Multiple-Range Test .......................-..+++- E-58 
17. SPSSX command file for Two-way ANOVA with repeated measurements .............. E-60 
18. SPSSX results for Two-way ANOVA with repeated measurements .....................-.... E-61 
19. SPSSX command file for Two-way ANOVA with more than one 

lace lita ace eet eeateatnemmarannenatnnnieniti E-64 
20. SPSSX results for Two-way ANOVA with more than one observation 

iia aan ita aiin ita anetadenearina ernment E-65 
21. STPK results for Two-way ANOVA with more than one observation 

STITT tiene ieiiiteeesileaiesidaa iad cen innntnemndnenseenantnenietiantatiabnanennetacenaseentnnenedl E-66 
ZZ. FPSB Command Tile Gor Cis-Saase 00Gl q....-.-ccccceecccecceecsecseceecceverssccncesscorecossocssoescoes E-70 
TTA TTT TNE E-71 
"TL ET Tt TEE T E-72 











Preface 


The information presented in this Technical Reference should be used in planning, analyzing, 
interpreting, and evaluating rangeland monitoring studies. 


‘There are important statistical considerations that must be borne in mind when planning and 
conducting monitoring studies of BLM rangelands. The recommendations resulting from 
these studies must be defensible both academically and legally. This requires a substantial 
background in the principles and procedures of statistical inference. 


Statistics is usually not the forte of the average BLM range conservationist, nor is it a subject 
matter that can be easily learned. There exists the feeling that all that is required is a few 
statistical formulas or a simple cookbook approach to assign meaning to data that were 
collected without regard for the type of data analysis to be performed. This feeling has been 
fostered in the Bureau by the frequent reference to “standard statistical procedures.” There 
has to be some understanding the statistical methods, and usually this is not the case when 
data are collected more or less arbitrarily and the conclusions “validated” statistically. 


Unfortunately, an average BLM range conservationist cannot become an accomplished 
statistician overnight, any more than a statistician can become an expert range conservationist 
by spending a couple of days in a BLM district. Understanding statistics takes time and 
effort, and there is no way around it. Once the statistical principles underlying the planning 
and conduct of monitoring activities arc understood, the specifics of monitoring design and 
anaiysis can be easily determined; however, trying to implement specifics without knowing 
the fundamentals easily leads to false conclusions and decisions that cannot be defended. 


The material that follows has been prepared to help BLM range conservationists cope with 
the statistical problems present in rangeland monitoring. The material is divided into five 
sections. Section 1 highlights the statistical topics required to analyze monitoring data and 
gives appropriate references. Section 2 addresses ‘he underlying statistical issues of rangeland 
monitoring. Sections 3 and 4 deal with the specific methods used for trend and utilization 
studies. Section 5 shows examples of data analysis. 

















Statistical Considerations in Rangeland Monitoring 
1. General Statistical Background 


This section does not address statistics in detail. Instead, the statistical topics required to 
analyze rangeland monitoring data are mentioned and references are given for the study of 
the individual statistical concepts and techniques. There is a large selection of introductory 
books used by universities in their range conservation curriculum to which the BLM range 
conservationist can refer for clarification of specific concepts. 


Characterization of vegetation data does not differ from the characterization of any other kind 
of data. Descriptive measures of location and variation are requized. Measures of location 
are the mean, median, and mode. Measures of variation or dispersion are the range, the 
variance, and the standard deviation. 


Inferences from the d.ata collected are the main concern of statistical analysis. The purpose of 
statistical inference is to draw conclusions about populations. A population is a collection 
of things that have some common observable characteristics and about which we want 
information. We can have populations of trees, plants, the weights of plants, and the ground 
cover produced by particular plants. 


In most cases, we cannot examine each element of the population and must be content with 
investigating only a part of the whole population. The part investigated is called a sample. A 
sample is a portion or subset of a population that is used to represent the population from 
which it is drawn. There are two types of samples: statistical samples and judgement 
samples. A statistical sample is one that is selected by a specific method of random selec- 
tion. A judgement sample is one that is picked from the target population based on the 
subjective decision of an individual. Statistical samples can provide an objective measure- 
ment of the population characteristics; judgement samples cannot. 


Sampling is done to draw inferences about the population sampled. The precision of infer- 
ences based on random samples can be assessed by means of probability theory. 


Several probability distributions are used for describing vegetation characteristics. Most 
statistical analyses are based on the normal distribution. Other commonly used distributions 
used in the analysis of vegetation data are the Poisson, binomial and log normal. 


Statistical analyses rest on the connection between sampling and inference. A statistic is a 
number computed from a sample. On the basis of the sample statistics, we make inferences 
about population parameters, which are numbers describing the content of the population 
under investigation. The most important fact to be used in statistical inference is that the 
distribution of sample means can be approximated by a normal distribution, and that the mean 
of this sampling distribution is the mean of the population. This is called the “central limit 
theorem”. 


A sample statistic estimates the value of the population parameter. The precision of such 
an estimate is measured by a quantity called the standard error, and this quantity is used to 
construct a confidence interval that we can be reasonably confident surrounds the unknown 
value of the population parameter. Confidence intervals for the population mean are easily 
estimated using the central limit theorem. 
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To compare two or more populations, we use tests of significance. A t-test is used to com- 
pare two sample means. Analysis of Variance (ANOVA) is used for comparing more than 
two sample means. Tests of significance involve the formulation of what is called a null 
hypothesis, an alternative hypothesis, and a level of significance. The level of significance 
gives the probability of rejecting the null hypothesis, when such a hypothesis is true. If we 
falsely reject the null hypothesis, then we are said to make a Type I error. Sometimes it is 
possible to accept a null hypothesis falsely, rather than reject it. In doing so we make a Type 
Il error. Balancing the two types of errors is very important in rangeland monitoring. 


All the above topics should be reviewed prior to continuing with Chapters 2 and 3 of this 
Technical Reference. Every range conservationist in the BLM has probably had at least one 
course in statistics during his/her college days. The textbook you used for such course should 
be a good reference. We recommend that you spend as much time as possible reviewing your 
old textbook. If your textbook is no longer part of your personal library, the following 
references are recommended: 


Elementary statistics: Glantz (1987); Huntsberger and Billingley (1987); Mattson (1986). 


More advanced statistics: Snedecor and Cochran (1980); Sokal and Rohlf (1981); Steel and 
Torrie (1980); Zar (1984). 


Elementary sampling theory: Scheaffer, Mendenhall, and Ott (1986); Williams (1978). 
More advanced sampling theory: Cochran (1977). 
Experimental design: Box, Hunter, and Hunter (1978); Winer (1971). 


All of the statistical procedures discussed and recommended in this Technical Reference 
(except for Chi-square analysis) are parametric; that is, they estimate population parameters 
(for example, means and variances) and test hypotheses concerning them. There are certain 
underlying assumptions that must be met or at least approximated in order for parametric 
statistical techniques to be valid. Chief among these is an underlying normal distribution and 
equality of variances. 


There is another set of statistical procedures available that do not require assumptions regard- 
ing distribution and variance. (They do, however, still require random sampling.) These 
procedures are termed nonparametric, or distribution-free, statistics. All of the parametric 
techniques described in this technical reference have nonparametric analogues. Many 
statisticians believe that, in general, parametric techniques are more powerful than nonpara- 
metric techniques, as long as the assumptions required to use parametric statistics are reason- 
ably well met. If in doubt about these assumptions, nonparametric methods can be used. 
Except for Chi-square analysis, which is a nonparametric technique, this Technical Reference 
does not discuss nonparametric methods. A considerable body of literature is available on 
this subject. Most of the statistics textbooks listed above discuss nonparametric statistics. In 
addition to these, the following books are valuable: Conover (1980); Lehmann (1975); and 
Sprent (1989). 


























Statistical Considerations in Rangeland Monitoring 
2. Monitoring of BLM Rangelands 


The BLM has several Technical References that deal with rangeland monitoring. (U.S. 
Department of the Interior, Bureau of Land Management, 1984a and 1985b). Technical 
Reference 4400-1, Planning for Monitoring, addresses the issue of sampling and sampling 
designs. Technical Reference 4400-7, Analysis, Interpretation and Evaluation, provides 
guidelines for data analysis. Both technical references give some general statistical guide- 
lines. 


Additionally, there are several well-known books on vegetation analysis and measurement 
techniques: Pieper (1978); Bonham (1989); Mueller-Dombois and Ellenberg (1974); Causton 
(1988); Krebs (1989); and Greig-Smith (1983). These books can help when planning a 
monitoring program. 

At the planning stage of a monitoring program, the questions to be answered, the hypotheses 
to be tested, and the effects to be estimated should be stated clearly. All this should be 
examined systematically. Clearly, there is no cookbook approach that can be followeu, but 


2.1 What are the objectives? 


The objectives should be stated clearly and precisely. Management objectives may be 
general, but the objectives of a sampling survey should be as specific as possible, each 
objective being stated as a hypothesis to be tested, a confidence interval to be computed, and 
a decision to be made. For example, the management objective may be to increase the 
percent cover of wheat grass from 10 to 15 percent in the Cranky Caribou key area. The 
objectives of the survey may be to estimate the percent cover with a margin of error of 2%, a 
90% confidence interval, and within certain personnel and budget constraints. 


2.2 What is sampled? 


The target population from which the sample is to be drawn and for which inferences are to 
be made must be defined. To define the target population, one must specify: 


a. The individual elements that make up the population. 
b. The attribute of interest. 
c. The way in which the attribute will be measured. 


As an example, suppose that one wants to estimate the percent cover on the Cranky Caribou 
key area. Several questions may be asked: 


a. Should any of the following areas be excluded? 
- fiparian area 
- areas with slopes of 3% or more 
- areas within 50 yards of established roads 
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b. What sort of cover is desired? 
- basal 
- canopy 
- ground cover, including rocks and litter 


c. How will the cover be measured? 
- line-intercept method 
- Daubenmire method 


Using the above example, the target population might be the basal area of all plants on the 
Cranky Caribou key area, as determined by the Daubenmire method on slopes of less 
than 3%. 


Next, the population must be divided into distinct sampling units which together constitute 
the population. Defining the sampling unit will answer the following questions: 


a. How will the individual units of the population be grouped together for sampling 
purposes? 
- individually (i.e., not grouped) 
- by plots or transects 
- by groups of plots or transects 


b. If by plots, what size and shape will be used? 
c. If by transects, how long will they be? 


Usually, the target population is individual plants, while the sampling unit is groups of plants, 
as defined by either transects or plots. Plots may be obtained by superimposing grids on 
maps of fields, forests, or other land areas. 


2.3 How is the sample drawn? 


The way in which samples are drawn is termed the sample design. There are several sample 
designs which may be appropriate for monitoring: simple random sampling, stratified 
sampling, and cluster sampling. The particular sampling design to be chosen is dictated by 
the nature of the problem and the availability of funds. 


The recommended approach for selecting monitoring sites is to randomly select sites in areas 
with similar management activities. Randomly allocating sampling units eliminates any 
possible bias on the part of the range conservationist, thereby increasing the accuracy of the 
analysis. The purpose of randomization is to guarantee the validity of the test of significance, 
based on the estimate of error made possible by repeated sampling. The number of samples 
and the sample design required for a given monitoring project depend mostly on the degree of 
precision desired. For a given sample design, the smaller the effect to b- detected, the greater 
the number of samples needed. 


Some consideration should be given to the question of interspersion when selecting monitor- 
ing sites. Interspersion refers to the distribution in space of the sampling units aiid is impor- 
tant for minimizing bias and the possibility of spurious effects. Randomization alone does 
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guarantee some degree of interspersion, but it may not be enough. If the randomized iayout 
appears segregated, the simplest and most widely used solution is to reject it and “re-random- 
ize” uniil a layout with an acceptable degree of interspersion is obtained. The importance of 
interspersion and its apparent contradiction to randomization are described in detail in 
Hurlbert (1984). As an aside, Hurlbert’s paper could be interpreted as justification for 
systematic sampling, at least in the sense that such sampling helps insure that the quadrants, 
plots, or plants measured are in interspaced throughout the sampling area. Prince (1986) 
offers some interesting insight on the issue of random versus systematic sampling of vegeta- 
tion, noting that systematic sampling .Aay be more accurate than a simple random sample of 
the same size because of the often strong correlation that exists between neighboring mea- 
surements in the natural environment. 


The method for data collection in the field is outlined in Technical Reference 4400-3 and 
Technical Reference 4400-4, (U.S. Department of the Interior, Bureau of Land Management, 
1984c and 1985a). As emphasized by Hunter (1980), high quality data can be obtained only 
if dependable measuring techniques are used and if trained, experienced, and reliable field 
personnel are employed to collect the data. 


For purposes of statistical analysis, it is preferable to randomly locate the plots or transects 
the first year and then measure the same plots or transects in succeeding years. Each plot or 
transect measured the first year is paired with the same plot or transect measured in subse- 
quent years. The data can be compared using a paired t-test for two years’ data and a re- 
peated measures analysis of variance for three or more years’ data. Since the members of 
each pair of plots or transects are positively correlated, tere is an increase in the ability of 
the statistical procedure to deiect a small difference. This is because the paired t-test and the 
repeated measures analysis of variance procedures eliminate a major source of variance, that 
existing from pair to pair. Instead of calculating the variance of differences among the 
individuals within each sample, the variance of the differences between the pairs is calcu- 
lated. Steel and Torrie (1980, pages 102-105) and the other statistics texts referenced in 
Section | elaborate on the advantages of using paired observations. 


2.4 Where are the samples taken? 


The field sample selection procedure is potentially the largest single source of (selection) 
bias. Samples selected because they are “convenient” or “representative” may save effort in 
some cases, but they may also produce seriously biased results. 


Random selection of sampling units can be easily accomplished. Let us illustrate this with an 
example. Suppose that we have determined that, for a given study, we need to select five 
transect locations in a key area. We proceed as follows: 


- Get a map of the key area. 

- Draw a grid directly on the map, or lay a grid over it. 

- Number all the grid squares from one to N. 

- Using a table of random num<is, or a random number generator on a computer or 
calculator, select five random numbers between one and N. This selects five 
locations on the map. 

- Check the spatial arrangement for interspersion. If not satisfactory, continue 
drawing sets of five random numbers until a satisfactory spatial arrangement is 
reached. 
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The five transect locations are now selected. Next, determine the method to locate the 
transect within the site. One could decide, for example, to proceed to the southwest corner of 
the site and walk northwest until a suitable transect location is reached. 


The size grid employed depends on: 


a. How accurately a spot on the n.ap can be located on the ground. There is no point in 
having a grid with 3,000 points if, in the field, only perhaps 300 points can be discerned. 


b. The number of potential sample (transect) sites. If there are only, say, 100 potential 
sites, it makes no sense to have nore than that number of grid points. 


Another method that can be employed to locate random coordinates in the field is described 
by Awbrey (1977). 


Yet another method of randomly locating transects is given in Technical Reference 4400-4 
(U.S. Department of the Interior, Bureau of Land Management, 19852), pages 31 and 32, and 
Illustration 24, page 101. This method involves permanently locating, a 100-foot baseline by 
means of two stakes. Transecis are then run perpendicular to the baseline at randomly 
selected points along a 100-foot tape. The direction of the transect (to the right or left of the 
baseline) is also determined randomly. Although the Technical Reference does not discuss it, 
care should be taken to run the baseline to minimize the amount of variability from one end to 
the other. If, for example, the area to be sampled is a hillside, anc the vegetation changes as 
one goes up or down the slope, the baseline should follow a contour line. The transects wili 
then cross the variability. The result will be that the transects will have the maximum 
amount of variability within them and the minimum amount of variability between them. 
This will make detection of differences more likely. 


What is probably not well understood concerning the use of the 100-foot baseline technique 
is that the target population consists only of those plants within an area that is 100 feet by 2 
times the length of each transect. For example, if each transect is 100 feet long, the target 
population is those plants within the 100-by-200-foot area that is potentially sampled by 
randomly placed transects. It is only for this area that statistical inferences can be made. Any 
inferences about the entire key area (assuming the key area is larger than 100 by 200 feet) are 


The number of points in a transect and the distance between points is an important consider- 
ation. Longer transects may show less variation than shorter ones, unless they cross vegeta- 
tion types, in which case they will show more variation than shorter ones. An increase in the 
size of the transect often results in a decrease in the number of replications that can be run, 
due to time and money constraints. Adequate replication of small transects is easier than 
adequate replication of large transects. It is usually more effi at to record more transects 
with fewer points per transect than fewer transects and more points per transect. 


The key areas sampled should be as homogeneous as possible if simple random sampling is 
used. Where key areas contain more than one vegetation type, stratified sampling should be 
used to reduve the total variance. If transects are the sampling units, they should cross the 
variability in the vegetation. Each sampling unit should have as much variability within it as 
possible and the variation between sampling units should be kept to a minimum. 
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2.5 How large a sample should be taken? 
The size of the sample required to detect differences depends on the following three items: 
a. The natural variability in the population. 
b. The precision required. 
c. The acceptable risk in the determination of confidence intervals. 


Large within-site variability requires more sampling units (unless one uses paired plots or 
transects—see discussion under 2.3 above). Often, some source of variation (such as eleva- 
tion or aspect) can be identified and the variation in the sampled variable can be at least 
partially explained. One can generally control the magnitude of this explained portion of the 
variation. The remaining unexplained variability will dictate the number of sampling units to 
take. 


The desired degree of precision indicates how close to the true mean the sample mean should 
be. Precision can be increased by increasing the sample size and by carefully selecting the 
sampling unit and sample design. 


The level of significance at which the null hypothesis is going to be tested determines the 
width of the confidence intervals. The level of significance is a specification of the accept- 
able risk that the actual confidence interval does not cover the true mean. The lower the risk 
tolerated, the larger the sample size. 


Setting the precision and confidence limits requires knowing the importance of the informa- 
tion being obtained and the resources available to do the necessary sampling and measure- 
ment. Some general guidelines are: 


a. Use 10% precision and 75% or 80% confidence when testing controversial hypoth- 
eses (i.c., when a future court appearance seems likely). 


b. Use 20% precision and 80% confidence for general resource management. 


The recommendation for a 75% or 80% level of significance versus the more conventional 
90% or 95% level is based on the need to lower the chances of committing a Type II error. 
Whenever we have controversy in the BLM we usually have resource deterioration. There- 
fore, we want to lower the chances of accepting the null hypothesis of no change in resource 
condition (Type Il error) when in fact the resource has deteriorated. The RISC (Range 
Inventory Standardization Committee) report recommends 80% as a reasonable confidence 
level. 


In summary, higher confidence, higher precision, and greater variability all dictate increased 
sample sizes. When an insufficient sample size is used, a significant difference may exist but 
not be apparent. Conclusions drawn from such an analysis may, therefore, be invalid. 


Technical References 4400-3 and 4400-4 (U.S. Department of the Interior, Bureau of Land 
Management, 1984c and 1985a), often recommend specific sample sizes. The sizes given in 




















Statistical Considerations in Rangeland Monitoring 


those documenis should be compared to the sample sizes determined to be necessary using 
the considerations discussed above for the particular level of significance, degree of preci- 

sion, and variability in the target population. Al! of these factors are combined in formulas 
for sample size determination (for example, Sokal and Rohif 1981, pages 262-264). 


An alternative to using formulas to determine sample size is to make use of sequential 
sampling techniques, whereby sampling continues until additiona! plots or transects do not 
significantly affect the mean (for the attribute of interest) of the more important (or abundant) 
plants. This technique is described by Mueller-Dombois and Ellenberg (1974, pages 77-80). 


2.6 What statistical tests should be used? 


Most range monitoring information is analyzed using Chi-square, confidence limits, t-test, 
Analysis of Variance (ANOVA), and time series regression analysis. The type of analysis 
selected depends upon (1) the number of plots or transects sampled, (2) the number of years 
of data, and (3) the type of data collected (cover, frequency, density). 


The two essential requirements of these tests are that (1) sampling must be random and (2) 
the sampling units (transects or plots) must be independent. Additional assumptions for 
ANOVA are (3) an underlying normal distribution and (4) equal variability within subgroups. 
Although plant attribute data and percent data usually violate requirement (3), data transfor- 
mation can be used to approximate a normal distribution and lack of normality can be easily 
mitigated. In general, lack of normality and inequality of subgroup variation tend to increase 
the overall variance, thus decreasing the chance of finding significant differences. 


If requirements (1) and (2) are violated, no valid statistical analysis is possible. Statistical 
analysis is not possible with data that have been arbitrarily collected. Statistical analysis is 
not just something added at the end of the investigative process to validate a foregone conclu- 
sion Or to dress up the results. On the other hand, if a sound sampling design has been used, 
the conclusions will have statistical validity and will be able to withstand outside scrutiny. 


Three common analyses of range monitoring data are (1) determining confidence limits for a 
given year, (2) comparing two years, and (3) comparing three or more years. 


1. Determining confidence limits for 2 species for a given year. 


- Statement you can make: “I am 90% confident that the percent cover of AGSP on 
the Cranky Caribou key area is between 13% and 21%” 
- Statistical test to use: Confidence limits on the year of interest. 


2. Comparing two years for a single species. 


- Statement you can make: “The increase in AGSP from 10% to 17% in the Cranky 
Caribou key area from 1983 to 1987 was statistically significant at the 15% level.” 

- Statistical test to use: t-test, paired if the same plots or transects were measured in 
the two years; two-sample if the same plots or transects were not remeasured. 
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3. Comparing three or more years for a single species. 


- Statement you can make: “There was a significant change in AGSP in the Cranky 
Caribou key area from 1975 to 1988,” or “there was a significant increasing linear 
trend in AGSP in the Cranky Caribou key area from 1975 to 1988.” 

- Statistical test(s) to use: If different plots were used each year, use analysis of 
variance and Duncan’s test or other multiple comparisons tests (Day and Quinn, 
1989) to determine which years are significanily different. The appropriate 
ANOVA model to use is that associated with a randomized complete block 
design, in which the treatments are the years and the blocks are the transects or 
plots (Winer, 1971). The analysis of variance can be used to detect linear and 
higher order trends. 


If the same plots weve re-measured in each year, the appropriate analysis of variance to use is 
that associated with a single-factor experiment having repeated measurer ents on the same 
elements. The treatments are the years, and a test for trend is used to detect a possible trend 
across the years (Winer, 1971). 


Once a statistical method is chosen to test the hypothesis, stick with the method. An unex- 
pected or undesired result is not by itself a valid reason for rejecting the method or the results. 
Examine all the steps involved in the analysis, beginning with the sampling design and the 
data collection phase. Scrutinize every aspect of the study. Then, make the decision whether 
to repeat the study or to use another sampling design and analysis. 


2.7 How are the statistical analyses run? 


All the statistical analyses required for rangeland monitoring can be run using existing 
computer packages, either on the [ioneywell mainframe or on a PC. 


The Honeywell mairframe in Denver has two general statistical packages and one specific 
rangeland monitoring program. The general statistical packages are STATPAK and SPSS. 
STATPAK can be accessed interactively and will suffice for t-tests, two-way analysis of 
variance, and linear regression. SPSS has to be accessed in a batch mode and it is a bit more 
cumbersome to use, but it is very comprehensive and is the best commercial package avail- 
able. The specific rangeland monitoring package is called A121/RENO/MONITOR (1985) 
and uses ANOVA and Duncan’s Multiple-Range Test to compare frequency and cover data 
across the years to detect any significant changes. The MONITOR program assumes a two- 
way randomized block design. The treatments are the years and the replications are the 
transects. The error term is the year x transect interaction. With the two-way randomized 
block design model, it is possible to identify significant differences between two or more 
years of data and two or more transects. In the instances when the number of transects is 
small and the two-way randomiz::d block design model shows no significant difference 
between years and between transects, a completely random design model with a single 
criterion of classification is also available to investigate possible significant differences 
between years. 
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There are many PC statistical packages. Just about any one will do mosi of the analyses 
required for rangeland monitoring. SPSS is available in the PC vzrsion and is called 
SPSSPC+. It comes in two packages: a basic package and an advanced package. They have 
the same capabilities as the corresponding mainframe SPSS packages. Other packages are 
STATGRAPHICS, MINITAB, SAS, ISP, SYSTAT, and others. 


2.8 Where and how should the information collected be stored? 


Information obtained from the monitoring program should be assembled in a format that is 
understandable by both resource specialist~ and decision makers. Document the objective(s) 
of the experiment, the details of the design of the experiment, the data collection phase, and 
the data analysis phase. 


The data collected should be entered into a data file or data base that is compatible with the 
statistical package used for the analysis. If the data are stored in a data base, updating and 
retrieval are considerably easier. Data base managers such as ASPEN/2 or dBASE III are 
suggested for data storage. ASPEN/2 is available on both the Honeywell mainframe and 

PCs, while dBASE III Plus is available on PCs. Since the data have to be entered into a 
computer for analysis, it makes sense to store the data in a data base format from which it can 
always be recalled for analysis or for further retrieval and update. 
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3. Analysis of Trend Study Methods 


The various trend study methods discussed in Technical Reference 4400-4 (U.S. Department 
of the Interior, Bureau of Land Management, 1985a) can be aggregated into three groups for 
the purposes of data preparation and data analysis. All of the methods discussed in the 
Technical Reference are covered here. If your method is slightly different from those dis- 
cussed in Technical Reference 4400-4, read the description of the methods covered by the 
three groups. Then select the closest group. 


3.1 Methods for photo plot, line intercept, and step point. 


This group includes the following methods from TR-4400-4: 














Data Collection Form 
Method Section Pages Illustration Pages 
Photo Plot 4.41 6-11 6 70-73 
Line Intercept 4.47 42-45 30 112-113 
Step Point 4.48 46-50 31 114-115 





Each plot or transect measured produces only one estimate of the percent cover, density, etc., 
for each species. The same plots or transects may or may not be remeasured in subsequent 
years. Plots may be of different sizes and transects may have different lengths from year to 
year. 


Obviously, if only one plot or transect per key area is established, the data obtained cannot be 
statistically analyzed. Several plots or transects (preferably 10 or more) are necessary for 
statistical analysis. In general, several shorter transects are better than one longer transect. 
Four 50-point transects will provide more information than one 200-point transect. 


The statistical analysis to be recommended is based on the following assumptions: (a) the 
plots or transects are randomly placed and (b) the transects are independent, while the plants 
or plots within transects are not independent. In addition, if the same plots are sampled every 
time, this has a bearing on the type of analysis used. 


3.1.1 Cover data preparation. 
a. Photo plot and line intercept methods. 


For each plot, use the percent cover computed on the data collection forms given in 
TR-4400-4. 


b. Step point method. 


Compute the percent cover for each transect for each species or cover category at each level 
(ground or basal and foliar, levels 1, 2, and 3). As an example, let us use the data on page 
115 of TR-4400-4 (U.S. Department of the Interior, Bureau of Land Management, 1985a) 
and work through the computations for ground and foliar I levels. This transect has only 100 
points rather than the 200 specified by the method, but the 100 points will make computations 
easier. 

11 














Statistical Considerations in Rangeland Monitoring 


All 100 points have to hit something on each level, even if that something is “thin air.” Thus 

the total hits for each cover layer must total 100. For the layers in which the count does not © 
add up to 100, one can assume that the missing values are the “thin air” hits. In the page 115 

example, 39 of the 100 points hit ground level without hitting anything above that layer, 

while 61 of the points hit some form of vegetation. To compute the percent cover at the 

ground level, combine the ground-level category cover with the basal or ground-level cat- 

egory from the basal and foliar cover as follows: 
































rumber percent 
category of hits cover 
B 16 16 
P 5 5 
N 23 23 
G 21 21 
C 6 6 
S 2 2 
R 0 0 
BOGR2 22 22 
SIHY 2 2 
BOCU 2 2 
OPUNT 1 1 
Total 100 100 
Compute the percent cover for the foliar - level 1 as follows: © 
number percent ~ 
category of hits cover 
BOCU 20 20 
BOGR2 3 3 
JUNIP 20 20 
ERIOG 2 2 
SIHY 2 2 
VIQUI 2 2 
N 4 4 
BOHI2 3 3 
P 2 2 
“Thin Air” 42 42 
Total 100 100 














Remember, there are 100 points and all of them must be accounted for. Therefore, those 
which are not listed on the data collection forms are the “thin air” hits. 
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3.1.2 Density data preparation. 
Compute the density by dividing the number of plants of each species by the square footage 


of the plot, i.c., 9 for a 3' x 3' plot and 25 fora 5'x 5' plot. Using the data on page 73 of TR- 
4400-4, the conversion to plants per square foot would be as follows for the listed species: 

















Number of Plants per 

Species plants scuare foot 
AGSP 4 4/9 = 0.44 
POSE *28 28/9 = 3.11 
SIHY 2 2/9 = 0.22 

| FIED 4 4/9 = 0.44 





3.1.3 Cover data file creation. ' 


Make separate data files for each species for each cover category you wish to compare. Note 
that cover category only applies to the Step Point Transect method, where more than one sort 
of cover data are collected. If only one cover category was collected, then only one data file 
for each species is necessary. Make separate columns for each year and separate rows for 
each transect. So, assuming that there are five transects, and that you wish to compare three 
years, the data file would have five rows and three columns. The data file for BOCU, for 
foliar - level 1, might look something like this: 








Transect 1984 1985 1986 
1 20 15 35 
2 24 14 6 
3 19 0 5 
4 5 13 14 
5 10 10 20 














3.1.4 Density data file creation. 


Make separate data files for each species you wish to compare. Make separate columns for 
each year and separate rows for each plot. The data file would have the same number of rows 
and columns as the one for percent cover, with the only difference being that the file would 
have density values rather than the cover percentages. 


3.1.5 Transformation of data. 


Transformation of cover data is required to make the data approximately normally distrib- 
uted. The recommended transformation is the arcsine square root, arcsinV y, with the cover y 
expressed as a decimal fraction. Transformation of the density data is recommended since 
enumeration data often follows a Poisson distribution. Use the square root transformation. 





1 Different statistical packages may require different data file structures than the one shown in this technical 
reference. 
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3.1.6 Analysis of one transect or plot per year. 
No statistical analyses are possible for either density and cover data. 
3.1.7 Analysis of two or more transects per year. 
a. Confidence intervals. 


To estimaie either the density or the percent cover for the key area in any given year, con- 
struct a confidence interval for the year of interest. 


b. Comparing two years. 


If the same plot locations were used in both years, use a paired t-test. If the same plot loca- 
tions were not used in both years, then use a two sample t-test. 


c. Comparing three or more years. 


If the same plot locations were used in all years, the appropriate analysis is an ANOVA for a 
single factor exneriment having repeated measurements on the same elements. If the same 
plot locations were not used in all years, ther an ANOVA for a randomized complete block 
design is the appropriate analysis. The treatments are the years and the blocks are the 
transects or plots. Duncan’s Multiple-Range Test can be used to detect individual differ- 
ences between the years. The ANOVA can also detect linear and higher-order trends be- 
tween years. 


3.2 Procedures for community structure analysis and the 
Daubenmire method. 


This group includes the following methods from TR-4400-4: 











Data Collection Form 
Procedure Section Pages Illustration Pages 
Community Structure 4.42 12-17 10-13 77-82 
Analysis 
Daubenmire Method 4.43 18-23 14-15 89-89 








Measurement of a plot or transect involves tallying the data into different groups so that the 
final figures are counts by category of the information obtained from reading the transect. 
The same plots or transects may or may not be measured in subsequent years. For analysis 
purposes, plots or transects can be of different sizes or lengths from year to year. 


The transects in these two methods should be permanently marked. However, for statistical 
analysis it will be assumed that new transects are established each year. This is because 
different people pace at different lengths, and the exact same spot of ground is not 
remeasured each year. Additional assumptions are that the transects are independent and 
randomly placed. The plots within the transects are usually not independent. 
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If only one transect per key area is established, no rigorous statistical analysis will be possible 
for either cover or frequency, since each transect produces only one data puint for cover or 
frequency. As pointed out earlier, several shorter transects make more sense than a longer 
transect. Four randomly placed, independent, 25-pace transects are more meaningful than a 
single 100-pace transect, because the replications make possible the estimation of sampling 
error. 


A questionable analysis for cover may be possible if only one transect per key area is estab- 
lished. A Chi-square analysis can be done to compare cover classes by years. The analysis 
will be suspect, however, because the data points along the transect were not collected in a 
truly random fashion and are probably not independent. Both of these conditions are required 
for a valid Chi-square test. 


When density data are being collected using the Community Structure Analysis and only one 
transect has been established, it may be possible to use statistical analysis. If the circular 
density plots within the transect are spaced far apart, their degree of interaction may be 
minimal. In that case, a 10-plot transect can be considered to have 10 independent density 
values. 


3.2.1 Cover data preparation. 


a. Ohne transect per year. 


The numbers of counts per cover class will be used as input to the statistical analysis. The 
appropriate analysis is the computation of the chi-square statistic from a two-way contin- 
gency table. This analysis requires at least five or more counts in each category, so combine 
cover classes to obtain five or more counts per category per year. As an example, let us use 
the data for 1984 on page 78, TR 4400-4, and make up the data for 1986 and 1988. Let us 
say that, for some reason, we are only interested in species beginning with the letter A, which 
will include ASTRA and ATOB. Also, let us combine all cover classes over 25% into one 
category. The table produced from these categories would look like this: 


Number of Counts by Cover classes for the “A” species. 








Cover class 1984 1986 1988 
1-5% 7 14 10 
5-25% 14 22 18 
25% and over 7 10 5 














Similar tables can be prepared for each species combination and cover class combination of 
interest. 


b. Two or more transects per year. 
If more than one transect per key area has been established, the data file construction and 
analysis will be quite different. One cover value is computed for each transect as described 


on page 85, TR-4400-4. These cover percentages will be used as the basis for the analysis. It 
makes no difference whether cover was computed using the 6- or 10-cover class method. 
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3.2.2 Density data preparation. 


First, convert the density count data to density per square foot (or sguare meter) by dividing 
each count by the size of the hoop used. For example, if a 9.6 squarc foot hoop is used, 
divide each count by 9.6 Converting the data allows comparison of data collected with 
different size hoops or plots. If species are combined, it makes no difference whether the 
species are added together before or after dividing by the area of the hoop or rectangle. 


For each species, or combination of species, of interest, make a data file containing the 
density per unit area values. Make columns for each year and rows for each plot. As an 
example, the data for HUA would be converted to density per square foot as follows: 











Plot count density=count/9.6 
1 4 42 
2 3 31 
3 0 .00 
4 7 .73 
5 2 21 
6 5 56 
7 7 73 
8 0 .00 
9 2 21 

10 6 63 











The zero counts must be included. However, because the cell counts are too small, this 
species should be combined with other species until the cell counts are greater than or equal 
to 5. Do combine species, but do not combine the plots since each plot is an independent 
estimate of the total. 


3.2.3 Frequency data preparation. 


If only one transect per year was established, no statistical analysis is possible. When several 
transects per year are run, the frequency data must be collected with the same size plot each 
year in order to obtain valid comparisons between years. 


3.2.4 Cover data file creation. 

a. One transect per year. 
Make separate data files for each combination of species and cover classes to be compared. 
Make separate columns for each year, and separate rows for each cover class or combined 
cover classes. The same cover classes must be used for all years being compared. 

b. Two or more transects per year. 
Make separate data files for each species or combination of species to be compared. Make 
separate columns for each year and separate rows for each transect. So, assuming that there 


are five transects and that one wishes to compare three years, the data file would have five 
rows and three columns. It might look something like tiis: 
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Transect 1984 1986 1988 
1 33 40 55 
2 24 54 66 
3 30 50 25 
4 55 33 44 
5 10 30 20 














The values in the table are the percent cover. 
3.2.5 Density data file creation. 


Make separate data files for each species or combination of species to be compared. If the 
density values computed for HLJA under the density data preparation section (TR-4400-4) are 
used for 1984, and contrived values are generated for 1986 and 1988, the data file could look 
like the following: 








Plot 1984 1986 1988 
l 42 40 55 
2 31 54 66 
3 .00 50 .25 
4 73 33 A4 
5 21 .30 .20 
6 56 66 45 
7 73 A3 16 
8 .00 22 43 
9 21 00 78 

10 63 21 00 














3.2.6 Frequency data file creation. 


As stated earlier, unless more than one transect per key area was measured, no statistical 
analysis is possible. If more than one transect was measured, the data file would be con- 
structed with one column for each year and one row for each transect, the same as for cover 
with two or more transects per year. Each species must have a separate data file. If species 
are combined, new percent freyuencies must be computed by combining the species on a 
plot-by-plot basis and then computing a frequency from those values. 


3.2.7 Transformation of data. 
No transformation of cover data is required if data on only one transect per year are available. 


Otherwise, cover data should be transformed using the arcsine square root transformation. 
Before applying the transformation, express the cover data as a decimal fraction. 
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Density data should be transfoi med using the square root transformation. Frequency data 
should be handled similarly to cover data. Use the arcsine square root transformation, 
expressing the frequency as a decimal fraction. 


3.2.8 Analysis of one transect per year. 


Two or more years of cover data car ". compared using a two-way contingency table and the 


Two years of density data can be compared using a t-test. Use a paired t-test if the samc plots 
were measured. Ot) rwise, use a two-sample t-test. 


If three or more years of density data are available, and if the same plots were remeasured, 
use ANOVA with repeated measurements. If the same plots were not remeasured, use a two- 
way analysis of variance in which the treatments are the years and the blocks are the plots. 
Duncan's Multiple-Range Test can be used to determine which years are different. 


No frequency comparisors are possible with only one transect per year. 
3.2.9 Analysis of two or more transects per year. 
a. Confidence intervals. 


To estimate the cover, density, or frequency for the key area in any ziven year, construct a 
confidence interval for the year of interest. 


b. Comparing two years. 


Cover and frequency data can be compared using a t-test. Use a paired t-test if the same plots 
were measured. Otherwise, use a two-sample t-test. 


To compare density data, either use a t-test, if comparing transects, or an ANOVA, if the 
transects are subdivided into plots. The analysis of variance model is a randomized block 
design with more than one observation per block. The treatments are the years. The blocks 
are the transects. The repeated observations are th= plots within the transects. 


c. Comparing three or more years. 


Cover and frequency data are analyzed using a two-way ANOVA for a randomized block 
design. The years are the treatments and the transects are the blocks. Duncan’s Multiple- 
Range Test can be used to detect differences between the years. The ANOVA of variance 
can also detect linear and higher-order trends. 


Density data can be analyzed two ways. If the plots within the transects are pooled, use the 
model recommended for cover and frequency. If the plots within the transects are treated as 
replicated measurements, use an ANOVA randomized block design with more than one 
observation per block. The treatments are the years, the blocks are the transects, and the 
repeated observations are the plots within the transects. 
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3.3 Methods for pace frequency, quadrat frequency, and nested 
frequency. 


This group includes the following methods from TR-4400-4: 














Data Collection Form 
Method Section Pages Illustration Pages 
Pace Frequency 4.44 24-28 18 93-94 
Quadrat Frequency 4.45 29-35 12 97-98 
Nested Frequency 4.46 36-41 25-28 102-109 





The sampling unit in the three methods is the transect. The transects are randomly placed and 
are assumed to be independent. For the purpose of the statistical analysis, it is assumed that 
new transects are established each year, since it is practically impossible to remeasure the 
same transect and plots year after year. The plots within the transects are not independent. 
Precision can be increased by increasing the number of transects. Accuracy is increased by 
increasing the number of quadrants within transects. 


3.3.1 Cover data preparation. 


Compute the percent cover for each transect for each cover category, similar to what is done 
for the total of the four transects on page 94, TR-4400-4. For example, here are the conver- 
sions for the cover categories for transect 1, page 94 of TR 4400-4: 

















Cover category count % cover 
Bare Ground 7 (7/50)x100= 14 
Persist. Litter 2 (2/50)x100= 4 
Non-Per. Litter 22 (22/50)x100= 44 
Rock 4 (9/50)x100= 8 
Live Veg. 10 (10/50)x100= 20 





Similar computations would be done for each transect for each cover category and each 
species. 


3.3.2 Frequency data preparation. 
Compute the frequency for eac)) plot size for each transect for each species, similar to what is 


done for the total of the four \:ansects on page 94 of TR-4400-4. For example, here are the 
conversions for the specie, | or transect 1, page 94 of the reference. 
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Species count frequency 
AGSP 25 (25/S0jx100 = SO 
KOCR 21 (21/S0)x100 = 42 
PONE3 40 (40/50)x100 = 80 
STV 16 (10/50)x100 = 
PHLO2 5 (5/S0)x100 = 
ASTER 8 (8/S0)x100 = 
PENST I (1/S0)x100 = 2 
GUSA2 3 (3/5S0)x100 = 6 
CHNA2 0 (0/S0)xi00 = O 
ARTR2 7 (7/50)x100 = 14 





Similar computations would be done for each transect for each cover category, each plot size, 
and each species. 


3.3.3 Cover data file creation. 


Make a separate data file for each cover category for which you wish to compare years. 
Make separate columns for each year and separate rows for each transect. So, assuming that 
there are four transects and that you wish to compar three years, the data file would have 
four rows and three columns. If the 1984 data for bare ground from page 94, in TR-4400-4, 
are used and we make up data for 1986 and 1988, the data file would look something like 
this: 

















Transect 1984 1986 1988 
l 14 04 25 
2 26 35 43 
3 52 39 65 
4 44 16 30 





3.3.4 Frequency data file creation. 


For nested frequency data, choose for each species the plot size that will keep the percent 
frequencies between 20 and 80 percent, if possible. Make a separate data file for each species 
and each plot size for which you wish to compare years. It will rarely be necessary to com- 
pare more than one plot size category for each species. For quadrat frequency, only one plot 
size is used, so a choice will not have to be made. 


Make separate columns for each year and separate rows for each transect. Assuming that 
there are four transects and that you wish to compare three years, the data file would have 
four rows and three columns. If the 1984 data for AGSP for plot size 3 from page 105, in 
TR-4400-4, are used for 1984 and data are made up for 1986 and 1988, the data file would 
look something like this: 


























Transect 1984 1986 1988 
| 36 50 52 
2 34 35 43 
3 40 39 65 
4 52 16 3 





Transformation of both frequency and cover data is required to make the data approximately 
normally distributed. The recommended transformation is the arcsine square root, arcsin y, 
with y expressed as a decimal fraction. 
3.3.6 Analysis of two or more transects per year. 

a. Confidence intervals. 
For both cover and frequency, do confidence intervals for the years of interest. 

b. Comparing two years. 


For both cover and frequency, use a two-sample t-test. The number of plots may be different 
for the two years being compared. 


c. Comparing three or more years. 
Cover and frequency data are analyzed using a two-way analysis of variance model for a 


Duncan's Multiple-Range Test can be used to detect differences between the years. The 
ANOVA can also detect linear and higher-order trends. 
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4. Analysis of Utilization Study Methods 


This section addresses all of the utilizatior, study methods described in Technical Reference 
4400-3. 


All utilizztion methods express utilization as a percentage of the current year’s production 
that has been removed. Each plot or transect measured produces an estimate of the percent 
utilization for each species, and several plots or transects are required to estimate the sam- 
pling error. 


The sampling unit in the paired plot method is the pair of plots. The difference between the 
protected and the unprotected plot represents the amount of forage consumed or destroyed 
during the foraging period. The protected and unprotected plots are dependent, but the pairs 
of plots are independent. 


In all of the other methods, the sampling unit is the transect. The transects must be indepen- 
dent and randomly placed. 


None of the utilization methods recommends that the same plots or transects be remeasured 
in subsequent years. 


As expressed repeatedly in the range monitoring section, there can be no statistical analysis if 
only one plot or transect is measured per year. No estimate of the sampling error is possible 
unless several samples are collected. Hence, several shorter transects are more meaningful 
than one longer transect. A 50-point transect will provide more information if it is separated 
into five transects of 10 observations each. 


4.1 Paired plot method data preparation. 


Compute the percent utilization for each plot for each species. The data for AGSP from TR- 
4400-3, illusization 2, page 54, will be used as an example. Use the same formula as for the 
overall percent utilization, (P-U)*100/P, but compute percent utilization for each plot as 
follows: 


PLOT! P-Ux 100 = 25-15x 100 = 10x 100 = 40% 
P 25 25 


PLOT 2 P-Ux 100 = 40-25x 100 = 15x 100 = 38% 
P 40 40 


PLOT3 P-Ux 100 = 30-15% 100 = 234.100 = 61% 
P 38 38 


PLOT4 P-Ux 100 = 30-18x 100 = 12x 100 = 40% 
P w 30 


PLOTS P-Ux 100 = 28-12% 100 = 16x 100 = 57% 
P 28 28 
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4.2 Data preparation in all the other methods. 


Compute the percent utilization for each transect and for each species as it is done in 
TR-4400-3. 


4.3 Data file creation for all methods. ' 


Make separate data files for each species you wish to compare. Make separate columns for 
each year and separate rows for each plot. Assuming that there are five plots or transects, and 
that you wish to compare three years, the data file would have five rows and three columns. 
It should look something like this: 








Transect 
or Plot 1984 1986 1988 
1 33 40 55 
2 24 54 66 
3 30 50 25 
4 55 33 44 
5 10 30 20 














4.4 Transformation of data. 


Use the arcsine square root transformation for all utilization data. Express utilizations as a 
decimal fraction. 


4.5 Analysis. 
a. Confidence intervals. 


Do confidence intervals on the years of interest. 


b. Comparing two years. 


Use a two sample t-test. This is not a paired t-test. If the same plot locations were used in 
both 1984 and 1988, then it would be a paired t-test. However, this should not be done. See 
Technical Reference 4400-3, page 7, (b) Continuing Study, for a discussion of this. Note also 
that, since it is a two-sample t-test, the number of plots may be different for the two years 


being compared. 
c. Compare three or more years. 


Use a two-way ANOVA for a randomized block design. The years are the treatments and the 
transects or plots are the blocks. Duncan’s Multiple- Range Test can be used to determine 
which years are different. The ANOVA can also detect linear and higher-order trends. 





2 Different statistical packages may require different data file structures than the one shown in this technical 
reference. 
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5. Examples of Data Analysis 


What follows are example problems that are solved using SPSSX and STPK on the 
Honeywell mainframe. SPSSX is used in a batch mode. STPK is used interactively. Each 
example is solved using both SPSSX and STPK. For SPSSX, the command file with the list 
of SPSSX commands required to solve the problem is shown first. The computer solution 
follows. For STPK, the whole interactive session is shown. 


The reason for showing both SPSSX and STPK is to show the reader that STPK, although a 
modest package compared to SPSSX, can accomplish most of the analyses required for 
rangeland monitoring. It does not have the capabilities to do Duncan’s Multiple-Range Test 
or linear contrasts, but it can do all of the other tests. Duncan’s can be done using A121/ 
RENO/MONITOR and linear (and higher-order) contrasts can be easily done by hand, as will 
be shown. 

SPSSX is a package that is geared to the statistician. Some may find it a bit too complex and 
may opt to use STPK or any of the other statistical packages available for the micro. By 


using STPK as a reference, you should be able to establish the suitability of other interactive 
statistical packages that you may have available on your PC. 


5.1 Data sets used. 


The following data sets are used in the examples that follow. 


ONECDATA. Cover count for 3 cover classes and 3 years. Cover classes are rows, years are 
columns. Exhibit 1. 


CHDATA. Cell entry counts by cover classes 1, 2, and 3, for years 84, 85, and 86. Exhibit 2. 


XXDATA. Frequency counts for 10 transects for 5 years. Cell entries by year and transect. 
Exhibit 3. 


YYDATA. Frequency counts for 10 transects for 5 years. Transects are rows, years are 
columns. Exhibit 4. 


DDDATA. Density for 2 transects, 4 years, and 3 replications. Cell entries by transect, year, 
and replication. Exhibit 5. 


5.2 Example of t-test. 
SPSSX. 
a. Independent Samples. 
The SPSSX command file is shown in Exhibit 6. This command file selects the data from 


data set XXDATA, transforms the data using the arcsin transformation, and compares years 
1 and 2 using a t-test for independent samples. The data are first divided by 20, since there 
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are 20 quadrants per transect, to get decimal fractions. Then, the arcsinv of the fraction is 
computed to obtain radians. The radians are converted to degrees by multiplying by 
360/n = 57.29578. 


The SPSSX results are shown in Exhibit 7. 
b. Paired Samples. 


The SPSSX command file is shown in Exhibit 8. This command file selects the data from 
data set YYDATA, transforms the data using the arcsinV transformation, and makes paired 
comparisons of year | and year 2 using a t-test for paired samples. The data transformation is 
the same as described in item a above. 


The SPSSX results are shown in Exhibit 9. 
STPK. 


The interactive session is shown in Exhibit 10. It uses data set YYDATA, transforms the 
data using the arcsin transformation. and compares year | and year 2 for both independent 
samples and paired samples using the t-statistic. The results agree with the SPSSX results. 


5.3 Example of Two-way ANOVA (Randomized block design). 


SPSSX. 


The SPSSX command file is shown in Exhibit 11. The command file uses data set 
XXDATA, transforms the data using the arcsinV transformation, and performs a Two-way 
ANOVA using a randomized block design. The years are the treatments and the transects are 
the blocks. It also computes linear and higher-order contrasts to detect trends. 


The SPSSX results are shown in Exhibit 12. The F value for YEAR indicates significant 
difference between the years at the 0.000 level of significance. The F value for the first 
PARAMETER under YEAR shows the significance of the linear trend at the 0.000 level of 
significance. Notice that the F values for the 2nd, 3rd, and 4th PARAMETER under YEAR 
(quadratic, cubic, and quartic trend components) are not significant at the 0.1 level. 


STPK. 


The interactive session is shown in Exhibit 13. It uses data set XXDATA, transforms the 
data using the arcsinV transformation, and performs a Two-way ANOVA. 


The linear and higher-order contrasts can be easily calculated by hand. We need to know the 
yearly totals for the transformed data and the contrast coefficients. The totals can be obtained 
by using the GROUP 1, TRANSFORMATION analysis option in the STPK menu, printing 
the data, and summing the totals for the years. The contrasts coefficients can be obtained 
from Winer (1971, p.878). 


For the XXDATA, the transformed year totals are: 
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Year 1 2 3 4 5 © 


444.00 537.16 562.50 652.46 767.64 

















and the contrast coefficients, c, are: 














Yc2 
Linear -2 -1 0 l 21 0 
Quadratic 2 -1 -2 -1 2 14 
Cubic -1 2 0 -2 1 10 
Quartic l 4 6 4 l 70 





Then, the linear contrast F can be computed (Winer, 1971, p.296): 


Chin = (-2)(444.00)-(1)(537.16)+(1)(652.46)+(2)(767.64) 
= 762.58 


Cyin = 581528.26 
MSjin = Cijn/n Dc” = 581528.26/((10)(10)) = 5815.2826 


Fin = MStin/MSresidual = 5815.2826/109.84 


= 52.94 } 


which agrees with the F for the linear trend in the SPSSX printout. 
Similarly for the quadratic F, 


Cquad = (2)(444.00)-(1)(537.16)-(2)(562.50)-(1)652.46)-(2)(767.64) 
= 108.66 


Cquad = 11806.996 

MS quad = Cquad/NXc? = 11806.996/((10)(14)) = 84.3356 
Fquad = MSquad/MSresidual = 84.3356/109.84 

= 0.7676 


which agrees with the F for the quadratic trend in the SPSSX printout. 
5.4 Example of Duncan’s Multiple-Range Test. 


SPSSX. 


The SPSSX command file is shown in Exhibit 14. This command file uses the data set @ 
XXDATA, transforms the data using the arcsin transformation, compares years using a One- 
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way ANOVA model, and conducts a Duncan’s Multiple-Range Test. The only reason for 
using the One-way ANOVA model is that the Duncan’s test has to be run with this model in 
SPSSX. The Duncan’s test result generated by this model should be good enough for our 
purposes, unless the mean square due to the transects is significant, in which case a manual 
computation can be easily performed to arrive at the Duncan’s statistic for the two-way 
model. 


The SPSSX results are shown in Exhibit 15. The results of Duncan’s test are best understood 
by comparing columns versus rows in the printout. Pairwise group comparisons of | vs. 2, 2 
vs. 3, and 3 vs. 4 are not significantly different, but 4 vs. 5 is. Pairwise comparisons of 1 vs. 
3, 2 vs. 4, and 3 vs. 5 are significantly different. Pairwise comparisons of 1 vs. 4, and 2 vs. 5 
are significantly different. Pairwise comparison of 1 vs. 5 is also significantly different. 


As indicated earlier, this should be good enough if the transect effect is not significant and, 
hence, the mean square error is pretty much the same for the One-way ANOVA and the Two- 
way ANOVA (in this example MS,,,o, = 116.63 for one-way and MS,.,5, = 109.84 for two- 
way). If the transect effect is significant, the significant differences in the Duncan’s test for a 
Two-way ANOVA model can be recomputed as follows: 
The SPSS printout gives this expression for the significant difference, S, 

S = 7.6364 x RANGE x V (1/N(i)+1/NGj)) 
which is equivalent to 

S = V (MSerroy/2) X RANGE x ¥ (2/N) 


Since the MSerror = 109.84 for the two-way model and N (number of transects) is 10, then, 
the significant difference expression becomes 


S = V (109.84/2) x RANGE x V (2/10) 


Using the RANGE values given in the printout, the values for S become 





2 3 4 5 











S 9.4455 9.9426 10.2409 10.4729 





Now, let us compare the means. The means are ordered in ascending order, 


l 44.40 
2 53.71 
3 56.25 
4 65.24 
5 76.76 
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For contiguous means, use Sz = 9.4455 
44.40 - $3.71 = 9.31 , not significant 
53.71 - 56.25 = 2.54 , not significant 
56.25 - 65.24 = 8.99 , not significant 
65.24 - 76.76 = 11.52 , significant 

For means two places apart, use S; = 9.9426 
44.40 - 56.25 = 11.85 , significant 
53.71 - 65.24 = 11.53 , significant 
56.26 - 76.76 = 20.51 , significant 

For means three places apart, use S, = 10.2409 


44.40 - 65.24 = 20.84 , significant 
53.71 - 76.76 = 23.05 , significant 


For means four places apart, use S; = 10.4729 
44.40 - 76.70 = 22.36 , significant 


Expressing this in the same graphical fashion as the SPSSX printout, 














Means 1 2 3 4 5 
l 
2 
3 
4 * 
5 * * * * 








Notice that the significant differences are the same as in the SPSSX printout, which is to be 
expected due to the lack of significance of the transect term, and the corresponding similar 
values of the mean square error terms in the one-way and the two-way models. 


STPK and A121/RENO/XMONITOR. 


STPK does not have the capability to run Duncan’s test. Al21/RENO/MONITOR does have 
the capability to conduct Duncan’s test on contiguous means as well as a Two-way ANOVA. 
The interactive session is shown in Exhibit 16, using data set YYDATA adapted to the format 
of the program. Notice that the mean squares, F values, and their level of significance for the 
years, transects, and residual error are identical to those obtained using either SPSSX or 
STPK for the Two-way ANOVA design. The values for S are slightly different, which 
appears to be due to the values used in building the Duncan tables in the computer subroutine. 
The “non-significant groupings” give the comparisons for contiguous means only. Notice 
that the results agree with the SPSSX printout. Pairwise comparisons of 84 vs. 85, 85 vs. 86, 
and 86 vs. 87 are not significantly different, but 87 vs. 88 is significantly different. 
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5.5 Example of Two-way ANOVA with repeated measurements. 


SPSSX. 


The SPSSX command file is shown in Exhibit 17. This command file uses data set 
YYDATA, transforms the data using the arcsin transformation, and performs a Two-way 
ANOVA with repeated measurements on the same plots or transects. Refer to Winer 
(1971, p.261) for a detailed explanation of the single-factor experimental design with re- 
peated measurements on the same element. 


The SPSSX results are shown in Exhibit 18. Although conceptually the repeated measure- 
ment design is different from the two-way randomized block design, the results are 
computationally the same. The within-cells mean square in the test of “between subjects” 
corresponds to the transect mean square. Similarly, in the test of “within subject,” the within- 
cells mean square corresponds to the residual mean square, and the year mean square is the 
same in both analyses. 


Since there is no advantage to using this design when dealing with only one factor (years) and 
the results are computationally the same, use a randomized block design even if the same 
plots or transects are remeasured. 


STPK. 


STPK does not have the capability to run this design. However, since the Two-way ANOVA 
design is computationally identical to the design of repeated measurements on the same plot 
or transect, use the two-way randomized block design. 


5.6 Example of Two-way ANOVA with more than one 
observation per block. 


SPSSX. 


The SPSSX command file is shown in Exhibit 19. This command file uses the data set 
DDDATA, transforms the data using the square root transformation, and performs a two- 
factor analysis of variance with repeated measurements on one factor. Refer to Winer 
(1971, p.518) for a detailed explanation of the two-factor experiment with repeated measure- 
ments on one factor. 


The SPSSX results are shown in Exhibit 20. The YEAR and the YEAR x TRANSect inter- 
action mean squares are tested against the residual mean square. The YEAR factor is signifi- 
cant at the 0.000 level, which means that there is a significant difference between the years at 
that level of significance. The YEAR x TRANSect interaction is only significant at the 0.166 
level. The TRANSect effect is tested against the ERROR | mean square. The F value shows 
significance at the 0.026 level. 
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STPK. 





The interactive session is shown in Exhibit 21. The analysis is conducted using a three-factor © 
design and pooling the appropriate mean squares to compute the two error terms. To exam- 

ine the significance of the year (Y) and the year x transect interaction (YT), we pool the YR 

and the T x YR term. This results ina POOL! value of 14.83, which is the same as the 

RESIDUAL in the SPSSX analysis. To examine the significance of the transect effect (T), 

we pool the replication (R) and the TR interaction. This results in a POOL2 value of 17.16, 

which is the same as the ERROR | term in the SPSSX printout. 


5.7 Example of Chi-square. 


SPSSX. 


The SPSSX command file is shown in Exhibit 22. This command file uses data set 
CHDATA, sets up a contingency table, and computes the Chi-square statistic. The hypoth- 
esis being tested is that of independence between the rows (cover classes) and the columns 
(years). This means that rows and columns represent an independent classification, that is, 
the proportion of each year total that belongs to each cover class is the same for all years. If 
this is the case, there is no difference between years and the Chi-square statistic will not be 
significant. If it is not the case, the years are different and the Chi-square statistic will be 
significant. 

The SPSSX results are shown in Exhibit 23. The Chi-square statistic obtained is only signifi- 
cant at the 0.8864 level. Since we are looking for significance at the 0.01 or 0.05 level, we © 
can state that the years are not significantly different at either of those two levels. 


STPK. 


The interactive session is shown in Exhibit 24. The analysis is conducted using the 
ONECDATA set. The Chi-square value obtained of 1.149 agrees with the SPSSX printout. 

















Key area. 


Hypothesis, Hp. 


Glossary of Terms 


The closeness of a measured or computed value to its true value. 


Probability of rejecting Ho, given that the Hp is true; that is, the 
probability of Type I error. 


A statement that there is a difference between the true value of the 
parameter of the population sampled and the hypothesized value of the 
population parameter. This is equivalent to saying that the treatment 
does have an effect. 


Probability of accepting Hg, given that the Ho is false; that is, the 
probability of Type II error. 


The process of making generalizations about a population from one or 
more samples. Statistical inference allows one to assess how reliable 
such generalizations are (in «ims of standard errors, confidence inter- 
vals, etc.). Logical inference does not allow an assessment of precision. 


A relatively small portion of a rangeland area selected because of its lo- 
cation, use, or grazing value and used as an area on which to monitor the 
effects of grazing use. It is assumed that key areas, if properly selected, 
will reflect the effects of current grazing management over all or a part of 
a pasture, allotment, or other grazing unit. Since key areas are subjec- 
tively located, statistical inferences concerning the population sampled 
are possible only for the key area (or the portion of the key area 
sampled). Any generalizations concerning the remainder of the pasture 
or allotmentare logical inferences. This does not mean that such gener- 
alizations are not valid, only that there is no way of assessing the degree 
of error associated with such generalizations. 


A statement that there is no difference (hence the term “null”) between 
the true value of the parameter of the population sampled and the 
hypothesized value of the population parameter. This is equivalent to 
saying that a treatment has no effect. 


A number computed from a population. 


Any set or collection of things that has some observable thing about 
which one wishes information. 

The closeness of repeated measurements of the same quantity. 

A portion or subset of a population. A collection of sampling units. 


Samples are used to represent the population from which the sample is 
drawn. 
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Sampling unit. A subdivision of the population that may be selected (through random 
sampling) through sampling. The sampling unit may be an individual 
plant, a transect, or a plot. The sampling units must cover the whole of 
the population to be sampled and not overlap. A number of sampling 
units selected through random sampling constitutes a sample. 


Statistic. A number computed from a sample. 


Type I error. The error of wrongly rejecting the null hypothesis, given that the null 
hypothesis is true (deciding that an ineffective treatment is effective). 

Type I error. The error of failing to reject the null hypothesis, given that the null 
hypothesis is false (failing to identify that an effective treatment is 
indeed effective). 


Type I and Type Il errors illustrated: 

















Ho 
Accepted Rejected 
Correct Type I 
True Decision Error 
Hoy 
Type I Correct 
False Error Decision 
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EXHIBIT 1. 


ONECDATA Data set 


7 14 10 
14 22 18 
7105 


EXHIBIT 2. 
CHDATA Data Set 


1 

2 

it 
8s 114 
85 2 22 
8S 3 10 
86 1 10 
86 2 18 


86 35 
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EXHIBIT 3. 


XXDATA Data Set 


1 
12 13 
13 13 
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46 16 
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49 16 
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7 12 10 16 20 
13 13 14 20 20 
13 10 12 16 14 
13 17 14 10 20 
8 8 17 11 18 

9 18 17 16 18 
7 141115 19 16 
8 7 14 10 18 17 
95 10 13 16 19 
10 9 15 15 17 20 
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EXHIBIT 4. 


YYDATA Data Set 


EXHIBIT 5. 


DDDATA Data Set 
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EXHIBIT 6. 





SPSSX Command File for t-test with independent samples 


10$$S,J7,T,MONI 

20$ : IDENT: SS01/FENNEC, SC210FKM 

30$ : SELECT: SPSS/X2.2/MAIN 

40$ : PRMFL:11,R,S,SC325FM/RANGEMON/XXDATA 
SOFILE HANDLE DATAINPT /UNIT=11 

60DATA LIST FREE FILE=DATAINPT 

70 /YEAR TRANS FREQ 

75SET WIDTH=80 

80LIST VARIABLES=ALL 

SOCOMPUTE FREQ=FREQ/20. 

1O00COMPUTE FREQ=ARSIN (SORT (FREQ) ) 
110COMPUTE FREQ=360.*FREQ/ (2.*3.1416) 
120T-TEST GROUPS=YEAR (1,2) / VARIABLES=FREQ/ 
180FINISH 

190$ : ENDJOB 
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EXHIBIT 7. © 


SPSSX results for t-test with independent samples 


GROUP 1 - YEAR EQ 1.00 
GROUP 2 - YEAR EQ 2.00 
VARIABLE NUMBER STANDARD STANDARD 
OF CASES MEAN DEVIATION ERROR 
FREQ 
GROUP 1 10 44.4009 9.375 2.965 
GROUP 2 10 53.7165 10.297 3.256 


* POOLED VARIANCE ESTIMATE SEPARATE VARIANCE ESTIMATE 


wy 

w * 

F 2-TAIL * T DEGREES OF 2-TAIL * T DEGREES OF 2-TAIL © 
VALUE PROB. * VALUE FREEDOM PROB. * VALUE FREEDOM PROB. 
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EXHIBIT 8. 





SPSSX command file for t-test with paired samples 


10$$S,J,T,MONI 

20$ : IDENT: SS01/FENNEC, SC210FKM 
30$:SELECT:SPSS/X2.2/MAIN 

40$ : PRMFL:11,R,S,SC325FM/RANGEMON/YYDATA 
SOFILE HANDLE DATAINPT /UNIT=11 
60DATA LIST FREE FILE=DATAINPT 

70 /TRANS FREQ] TO FREOQS 

75SET WIDTH=80 

8O0LIST VARIABLES=ALL 

85DO REPEAT FREQ=FREQ1 TO FREQS 
SOCOMPUTE FREQ=FREQ/20. 

1OOCOMPUTE FREQ=ARSIN (SORT (FREQ) ) 
110COMPUTE FREQ=360.*FREQ/(2.*3.1416) 
115END REPEAT 

120T-TEST PAIRS=FREQ1 FREQ2 
180FINISH 

190$ : ENDJOB 
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EXHIBIT 9. 


SPSSX results for t-test with paired samples 


DEGREES OF 2-TAIL 
FREEDOM PROB . 


ie + ee be. + 2 8 --T-TEST------ 
VARIABLE NUMBER STANDARD STANDARD 
OF CASES MEAN DEVIATION ERROR 
FREQ1 
10 44.4009 9.375 2.965 
10 53.7165 10.297 3.256 
FREQ2 
(DIFFERENCE) STANDARD STANDARD * 2-TAIL * T 
MEAN DEVIATION ERROR * CORR. PROB. * VALUE 
~9.3156 12.960 4.098 * 0.135 0.711 * -2.27 


4} 
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EXHIBIT 10. 





STPK results for t-tests with independent and paired samples 


TO 
TYPE A NUMBER FROM 1 TO 5S, SOS, OR 


TO CHANGE TO A DIFFERENT STAT GROUP, TYPE FINISH IN ANSWER TO 
-WHAT ANALYSIS DO YOU WISH TO PERFORM- 


TO TERMINATE STATPKG, TYPE END IN ANSWER TO 
-WHICH STAT GROUP WOULD YOU LIKE TO USE- 


IS YOUR DATA COMING FROM OTHER THAN TTY * 
=YES 

ENTER CATALOG FILE STRING FOR DATA FILE 
=SC325FM/RANGEMON/YYDATA 


DO YOU WISH TO LIST NAMES OF AVAILABLE ANALYSES * 
=YES 


THE FOLLOWING ANALYSES ARE AVAILABLE IN THIS GROUP.. 





BET. cece cccccccccccccccs ADD, REPLACE, AND DELETE DATA 
TRANSFORMATION........... LOG, SQUARE ROOT, ETC 
ELEMENTARY STATISTICS MEAN, STANDARD DEVIATION, ETC 
CORRELATION 

CROSS TABULATION......... FREQUENCY TABLE FOR 2 CROSSED VARIABLES 
SCATTER DIAGRAM 

HISTOGRAM 

LINE PLOT..........-0006% VERTICAL LINE PLOT 

RANK CORRELATION 

CHI-SQUARE... ......-55065 CONTINGENCY TABLE 

BD FEE occ ccccccccccccsces TWO SAMPLE MEANS 


WHAT ANALYSIS DO YOU WISH TO PERFORM * 
= TRANS FORMATION 


DO YOU WISH TO READ NEW DATA * 
=YES 


HOW MANY ROWS IN YOUR DATA MATRIX * 
#10 
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EXHIBIT 10 (CONT.) 


HOW MANY COLUMNS * 


=6 


SPECIFY INPUT CODE FOR YOUR DATA * 
=DF 


THE DATA WILL NOW BE READ FROM DISK IN A FREE FORMAT 


DO YOU WISH TO PRINT THE DATA JUST READ IN 


=YES 
1.000 7.000 12.000 10.000 16.000 20.000 
2.000 13.000 13.000 14.000 20.000 20.000 
3.000 13.090 10.000 12.000 16.000 14.000 
4.000 13.000 17.000 14.000 10.000 20.000 
5.000 8.000 8.000 17.000 11.000 18.000 
6.000 9.000 18.000 17.000 16.000 18.000 
7.000 14.000 11.000 15.000 19.000 16.000 
8.000 7.000 14.000 10.000 18.000 17.000 
9.000 5.000 10.000 13.000 16.000 19.000 
10.000 9.000 15.000 15.000 17.000 20.000 


DO YOU WISH TO CHANGE SOME VALUES * 
=NO 


TYPE TRANSFORMATION CODE * 
=10 


SPECIFY VARIABLE X AND CONSTANT, AND THE VARIABLE NUMBER ASSIGNED TO 
THE TRANSFORMED DATA IN THE FORM XX/XX.X/XX 
=2/0.05/7 


TYPE TRANSFORMATION CODE * 
=10 


SPECIFY VARIABLE X AND CONSTANT, AND THE VARIABLE NUMBER ASSIGNED TO 
THE TRANSFORMED DATA IN THE FORM XX/XX.X/XX 
=3/0.05/8 


TYPE TRANSFORMATION CODE * 
=1 


SPECIFY VARIABLE X AND THE VARIABLE NUMBER ASSIGNED TO THE TRANSFORMED 
DATA IN THE FORM XX/XX 
=7/9 


TYPE TRANSFORMATION CODE * 
=1 


“A 


EXHIBIT 10 (CONT.) 





SPECIFY VARIABLE X AND THE VARIABLE NUMBER ASSIGNED TO THE TRANSFORMED 
DATA IN THE FORM XX/XX 
=8/10 





TYPE TRANSFORMATION CODE * 
=12 


SPECIFY VARIABLE X AND THE VARIABLE NUMBER ASSIGNED TO THE TRANSFORMED 
DATA IN THE FORM XX/XX 
=9/11 


TYPE TRANSFORMATION CODE * 
=12 


SPECIFY VARIABLE X AND THE VARIABLE NUMBER ASSIGNED TO THE TRANSFORMED 
DATA IN THE FORM XX/XX 
=10/12 


TYPE TRANSFORMATION CODE * 
=10 


SPECIFY VARIABLE X AND CONSTANT, AND THE VARIABLE NUMBER ASSIGNED TO 
THE TRANSFORMED DATA IN THE FORM XX/XX.X/XX 
=11/57.29578/13 


TYPE TRANSFORMATION CODE * 
=10 





SPECIFY VARIABLE X AND CONSTANT, AND THE VARIABLE NUMBER ASSIGNED TO 
THE TRANSFORMED DATA IN THE FORM XX/XX.X/XX 
=12/57.29578/14 


TYPE TRANSFORMATION CODE * 
=0 


HOW MANY VARIABLES DO YOU WISH TO RETAIN * 
#14 


DO YOU WISH TO PRINT THE DATA MATRIX * 
=NO 


WHAT ANALYSIS DO YOU WISH TO PERFORM * 
=T TEST 


DO YOU WISH TO USE THE SAME DATA * 
=YES 
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EXHIBIT 10 (CONT.) 


SPECIFY CODE OF DESIRED HYPOTHESIS * 
=SOS 


THE FOLLOWING CODES SIGNIFY HYPOTHESES AVAILABLE FOR TESTING... 

1 MEAN OF A = GIVEN VALUE 

2 MEAN OF A = MEAN OF B, ASSUMING TWO VARIANCES ARE EQUAL 

3 MEAN OF A = MEAN OF B, ASSUMING TWO VARIANCES ARE NOT BQUAL AND 
SAMPLES ARE SMALL 

4 MEAN OF A = MEAN OF B, WHERE A AND B ARE RELATED 


SPECIFY CODE OF DESIRED HYPOTHESIS * 
=2 


SPECIFY THE COLUMN NUMBERS OF VARIABLES A AND B IN THE FORM AA/BB 


=13/14 
COMPUTED T VALUE...... 2.11535 
DEGREES OF FREEDOM... . 18 


DO YOU WISH TO USE THE SAME DATA BASE TO COMPUTE MORE T VALUES 
=YES 


SPECIFY CODE OF DESIRED HYPOTHESIS * 
=4 


SPECIFY THE COLUMN NUMBERS CF VARIABLES A AND B IN THE FORM AA/BB 


=13/14 
COMPUTED T VALUE...... 2.27302 
DEGREES OF FREEDOM... . 9 


DO YOU WISH TO USE THE SAME DATA BASE TO COMPUTE MORE T VALUES 
=NO 


WHAT ANALYSIS DO YOU WISH TO PERFORM * 
=FINISH 


WHICH STAT GROUP WOULD YOU LIKE TO USE, 
TYPE A NUMBER FROM 1 TO 5S, SOS, OR END. * 
=END 
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yd] Exhibit 11. 


SPSSX command file for two-way ANOVA with randomized block design 


10$$S,J,T,MONI 

20$ : IDENT: SS01/FENNEC, SC210FKM 

30$ : SELECT: SPSS/X2 .2/MAIN 

40$ : PRMFL:11,R,S, SC325FM/RANGEMON/XXDATA 
SOFILE HANDLE DATAINPT /UNIT=11 

60DATA LIST FREE FILE=DATAINPT 

70 /YEAR TRANS FREQ 

75SET WIDTH=80 

80LIST VARIABLES=ALL 

SOCOMPUTE FREQ=FREQ/20. 

1OOCOMPUTE FREQ=ARSIN (SORT (FREQ) ) 
11OCOMPUTE FREQ=360.*FREQ/ (7% .*3.1416) 
120MANOVA FREQ BY TRANS(1,10) YEAR(1,5) / 
130 PRINT=OMEANS (TABLES (CONSTANT YEAR) ) / 
140 PRINT=SIGNIF (SINGLEDF) / 

150 PRINT=PARAMETER (ESTIM) / 

160 CONTRAST (YEAR) =<POLYNOMIAL/ 

170 DESIGN CONSTANT YEAR TRANS/ 
180FINISH 

190$ : ENDJOB 
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EXHIBIT 12. © 


SPSSX results for two-way ANOVA with randomized block design 


e*eee2e2eCQARaALYysir Ss OP VARIANCE -- DEBSIGN i**** * * 
COMBINED OBSERVED GRAND MEANS 
VARIABLE .. FREQ 


WGT 59 .27571 
UNWGT . 59 .27571 





1 WGT. 44 .40092 
UNWGT . 44 .40092 
2 WGT. 53.71648 
UNWGT . 53.71648 
3 WoT. 56.25011 
UNWGT . 56.25011 
4 WGT. 65 .24654 
UNWGT . 65 .24654 
5 WGT . 76.76449 
UNWGT . 76 .76449 
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EXHIBIT 12 (CONT.) 


OF 


VARIANCE -- 


DESIGN 


TESTS OF SIGNIFICANCE FOR FREQ USING UNIQUE SUMS OF SQUARES 
SOURCE OF VARIATION 


RESIDUAL 


ESTIMATES FOR FREQ 


175680 


3954.31 
-48 
6028.27 
5815.16 
84 .36 
86.55 
42.19 
1294.05 
114.71 
233.38 
326.20 
10.33 
282.80 
46.51 
2.00 
114.47 
163.64 


CONSTANT 

PARAMETER COEFF. STD. ERR. 
1 59 .2757089 1.48217 

YEAR 
PARAMETER COEFF. STD. ERR. 
2 24 .1146448 3.31424 
3 2.90447908 3.31424 
4 2 . 94200420 3.31424 
5 -2.0541185 3.31424 

TRANS 
PARAME™ COEFF. STD. ERR. 
Ct) 1809171 4.44652 
7 9.57346855 4.44652 
8 -§.3315712 4.44652 
9 3 .27042337 4.44652 
10 -6.2535970 4.44652 
11 3 .90593692 4.44652 
12 1.75867973 4.44652 
13 -3.9080702 4.44652 
14 -5 .4272780 4.44652 


DF 
36 
1175 
4 i 
i Ss 
1 
1 
i 
9 
1 
i 
1 
1 
1 
1 
1 
1 
i 


T-VALUE 


39 .99242 


T-VALUE 


7.27607 
.87636 
.88769 

-.61979 


T- VALUE 


-.49048 
2.15302 
-1.19904 
. 73550 
-1.40640 
.87643 
.39552 
-. 87891 
-1.22057 


109. 
680. 
507. 
815. 
84. 
86. 
42. 
143. 
114. 
233. 
326. 
10. 
282. 
46. 
2. 
114. 
163. 


INDIVIDUAL UNIVARIATE .9500 CONFIDENCE INTERVALS 


F SIG OF F 

1599.39 .000 
13.72 .000 
52.94 .000 

-77 .387 

.79 -381 

.38 .539 

1.31 .266 

1.04 .314 

2.12 .154 

2.97 .C93 

.09 . 761 

2.57 .117 

-42 .519 

.02 .893 

1.04 .314 

1.49 .230 

IG. T LOWER -95% CL- 
00000 56.26972 62 


SIG. T LOWER -95% 


.00000 
. 38664 
. 38060 
.-53930 


17.39305 
-3.81711 
-3.77959 
-8.77571 
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UPPER 


.28170 


UPPER 


. 83624 
.62607 
-66360 
.66747 


UPPER 


.83705 
.59143 
.68639 
. 28839 
. 76437 
.92390 
. 77664 
.10989 
.59068 








EXHIBIT 13. 





STPK results for two-way ANOVA with randomized block design 


=3 


TO CHANGE TO A DIFFERENT STAT GROUP, TYPE FINISH IN ANSWER TO 
-WHAT ANALYSIS DO YOU WISH TO PERFORM- 


TO TERMINATE STATPKG, TYPE END IN ANSWER TO 
-WHICH STAT GROUP WOULD YOU LIKE TO USE- 


IS YOUR DATA COMING FROM OTHER THAN TTY * 

=YES 

ENTER CATALOG FILE STRING FOR DATA FILE 

=SC325 /RANGEMON /XXDATA 

UNABLE TO ACCESS FILE AS REQUESTED PLEASE RENTER 
ERROR STATUS = 400100000000 

ENTER CATALOG FILE STRING FOR DATA FILE 
=SC325FM/RANGEMON /XXDATA 


SPECIFY THE NUMBER OF FACTORS * © 
=2 


SPECIFY ONE-LETTER LABELS FOR EACH FACTOR IN THE FORM X/X/X/... 
=Y¥/T 


SPECIFY LEVELS (CATEGORIES) FOR EACH FACTOR IN THE FORM X/X/X/... 
=5/10 


SPECIFY INPUT CODE FOR YOUR DATA * 
=DF 


EACH INPUT RECORD SHOULD CONTAIN A SINGLE OBSERVATION PRECEDED BY ITS 
IDENTIFING SUBSCRIPTS - ONE FOR EACH FACTOR 


THE DATA WILL NOW BE READ FROM DISK IN A FREE FORMAT 


DO YOU WISH TO PRINT THE DATA JUST READ IN 
«YES 











EXHIBIT 13 (CONT.) 
FACTOR LEVELS OBSERVATION 
1 7.00000 
1 12 .00000 
1 10.00000 
1 16 .00000 
1 20 .00000 
2 13 .00000 
2 13 .00000 
2 14 .00000 
2 20 .00000 
2 20.00000 
3 13 .00000 
3 10.00000 
3 12 .00006 
3 16.00000 
3 14.00000 
4 13 .00000 
4 17.00000 
4 14.00000 
4 10.00000 
4 20.00000 
5 8.00000 
5 8.00000 
5 17.00000 
5 11.00000 
5 
6 
6 
6 
6 
6 
7 
7 
7 
7 
7 
8 
8 
8 
8 
8 
9 
9 
9 
9 
9 





18 .00000 
9.00000 
18 .00000 
17.00000 
16 .00000 
18 .00000 
14.00000 
11.00000 
15 .00000 
19 .00000 
16 .00000 
7.00000 
14.00000 
10 .00000 
18 .00000 
17.00000 
5.00000 
10.00000 
13.00000 
16 .00000 
19 .00000 
10 9.00000 
10 15 .00000 
10 15 .00000 
10 17.00000 
10 20.00000 
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EXHIBIT 13 (CONT.) 


DO YOU WISH TO MAKE ANY CORRECTIONS TO THE DATA * 


=NO 


DO YOU WISH TO PERFORM DATA TRANSFORMATION 
=YES 


SPECIFY TRANSFORMATION CODE * 
=10 


SPECIFY THE VALUE OF THE CONSTANT 
0.05 


SPECIFY TRANSFORMATION CODE * 
=1 


SPECIFY TRANSFORMATION CODE * 
=12 


SPECIFY TRANSFORMATION CODE * 
=10 


SPECIFY THE VALUE OF THE CONSTANT 


=57.29578 
SPECIFY TRANSFORMATION CODE * © 
=0 
FACTOR LEVEL 
Y 5 
T 10 
GRAND MEAN 59 .27585 
LINE SOURCE OF SUM OF DEGREES OF MEAN 
NUMBER VARIATION SQUARES FREEDOM SQUARE 
1 Y 6028 .29730 4 1507.07433 
2 T 1294 .05168 9 143.78352 
3 YT 3954 .32831 36 109 .84245 
TOTAL 11276 .67737 49 


DO YOU WISH TO POOL COMPONENTS * 
=NO 


DO YOU WISH TO COMPUTE F-VALUES * 
=YES 
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EXHIBIT 13 (CONT.) 





ERROR TERM IN THE FORM XX/XX. TYPE 0/0 WHEN YOU WISH TO STOP 
=1/3 


F-VALUE..... 13.72033 (D.F.= 4, 36) 


TYPE LINE NUMBERS FOR NEXT F CALCULATION 
=2/3 


TYPE LINE NUMBER OF MAIN EFFECT OR INTERACTION, AND LINE NUMBER OF 
F-VALUE..... 1.30900 (D.F.= 9, 36) 

TYPE LINE NUMBERS FOR NEXT F CALCULATION 

=0/0 

=NO 


DO YOU HAVE OTHER ANALYSIS-OF-VARIANCE PROBLEMS TO SOLVE 
=NO 


WHICH STAT GROUP WOULD YOU LIKE TO USE, 
TYPE A NUMBER FROM 1 TO 5, SOS, OR END. * 
=END 


DO YOU WISH TO POOL MORE COMPONENTS OR COMPUTE MORE F-VALUES 

















EXHIBIT 14. © 


SPSSX command file for Duncan’s multiple range test 


10$$S,J,T,MONI 

20$ : IDENT: SS01/FENNEC, SC210FKM 

30$ : SELECT: SPSS/X2.2/MAIN 

40$ : PRMFL:11,R,S,SC325FM/RANGEMON/XXDATA 
SOFILE HANDLE DATAINPT /UNIT=11 
60DATA LIST FREE FILE=DATAINPT 

70 /YEAR TRANS FREQ 

7SSET WIDTH=80 

S8OLIST VARIABLES=ALL 

SOCOMPUTE FREQ=FREQ/20. 

1LOOCOMPUTE FREQ=ARSIN (SORT (FREQ) ) 
L1OCOMPUTE FREQ=360.*FREQ/ (2.*3.1416) 
120ONEWAY FREQ BY YEAR(1,5) / 

150 POLYNOMIAL = 2/ 

160 RANGES = DUNCAN/ 

180FINISH 

190$ : ENDJOB 





EXHIBIT 15. 





SPSSX results for Duncan’s multiple range test 


cee ee ee eee ee eee eC OMBWAY- - - - - - ee ee ee eee 


ANALYSIS OF VARIANCE 


SUM OF MEAN F F 
SOURCE D.F. SQUARES SQUARES RATIO PROB. 
BETWEEN GROUPS 4 6028 .2688 1507 .0672 12.9218 .0000 
LINEAR TERM 1 5815 .1609 5815 .1609 49.8599 .0000 
DEVIATION FROM LINEAR 3 213.1079 71.0360 .6091 .6126 
QUAD. TERM 1 84.3600 84.3600 .7233 .3996 
DEVIATION FROM QUAD. 2 128.7479 64.3740 .5519 .5797 
@ WITHIN GROUPS 45 5248 .3554 116.6301 
TOTAL 49 11276 .6242 
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EXHIBIT 15 (CONT.) © 


“fa © oo ee 2 ----+--ONMBWAY--- +--+ ee - ee oe eg 


MULTIPLE RANGE TEST 


DUNCAN PROCEDURE 
RANGES FOR THE 0.050 LEVEL - 


2.85 3.00 3.09 3.16 
THE RANGES ABOVE ARE TABLE RANGES. 
THE VALUE ACTUALLY COMPARED WITH MEAN(J) -MEAN(I) IS.. 
7.6364 * RANGE * DSORT(1/N(I) + 1/N(J)) 


(*) DENOTES PAIRS OF GROUPS SIGNIFICANTLY DIFFERENT AT THE 0.050 LEVEL 


GGGGG 
RRRRR 
PPPPP 
MEAN GROUP 12345 © 
44.4009 GRP 1 
53.7165 GRP 2 
56.2501 GRP 3 ° 
65.2465 GRP 4 ee 
76.7645 GRP 5 secs 


HOMOGENEOUS SUBSETS (SUBSETS OF GROUPS, WHOSE HIGHEST AND LOWEST MEANS 
DO NOT DIFFER BY MORE THAN THE SHORTEST 
SIGNIFICANT RANGE FOR A SUBSET OF THAT SIZE) 


SUBSET 1 
GROUP GRP 1 GRP 2 

MEAN 44.4009 53.7165 

SUBSET 2 

GROUP GRP 2 GRP 3 ‘ 
MEAN 53.7165 56.2501 
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© EXHIBIT 15 (CONT.) 


SUBSET 3 

GROUP GRP 3 GRP 4 
MEAN 56.2501 65.2465 
SUBSET 4 

GROUP GRP 5 

MEAN 76.7645 
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EXHIBIT 16. 


/RENO/XMONITOR results for Duncan’s multiple range test 


SC325FM/RENO/XMONITOR 
INDICATE THE TYPE OF DATA TO BE ENTERED 
BY TYPING: 
1 - FOR FREQUENCY AND/OR COVER DATA 
2 - FOR DENSITY DATA 
=1 
DO YOU WANT TO UPDATE THE DATA BASE? (YES,NO OR NEW) 
=NO 
ENTER NAME OF DATA BASE FILE 
=SC325FM/RANGEMON/ZZBASE 
DO YOU WANT TO PRINT THE DATA BASE(YES OR NO)? 
=YES 


130.0 000000000000 0000000000000000000000000 Gael 
GAAAAAAXX 5 

Y¥8401011020 7131313 8 914 759 
¥850101102012131017 81811141015 
Y¥860101102010141214171715101315 
¥870101102016201610111619181617 
¥880101102020201420181816171920 


ENTER PERCENT LEVEL OF SIGNIFICANCE(S OR 10) 

25 

ENTER KEY AREA NUMBER; 

IF ALL KEY AREAS IN ALLOTMENT, ENTER ’ ALL’ 

=ALL 

ENTER ‘F’ FOR FREQUENCY ANALYSIS; ‘C’ FOR COVER. 

-F 

DO YOU WANT THE DETAILED STATISTICAL PRINTOUT? (YES OR NO) 
=YES 

IF YOU WANT THE OUTPUT TO GO TO A PERM FILE, ENTER FILE NAME; 
ELSE ENTER ‘NO’ 

=NO 











EXHIBIT 16 (CONT.) 





BLM ADMIN UNIT XXXXXXXX ;WILDLIFE XXXXX ;BIG GAME XXXXxXxXxxX 
ALLOTMENT XXXX ; PASTURE XX ;KEY XXxXxX 
PLANT SPECIES AAAAAA 


TWO-WAY ANOVA RESULTS 


YEAR RAW MEANS TRANS. MEANS 

84 49.00 44.4010 

85 64.00 53.7166 

86 68.50 56.2502 

87 79.50 65.2467 

88 91.00 76.7647 
SOURCE DP ss MS F PROB >F 
YEARS 4 6028 .2910 1507 .0728 13.72 0.0000 
TRANS 9 1294 .0473 143.7830 1.31 0.2665 
ERROR 36 3954 .3355 109 .8427 
TOTAL 49 11276 .6738 
STANDARD ERROR OF THE MEAN (TRANFORMED) 3.31 





DUNCANS MULTIPLE RANGE TEST AT THE 0.05 SIGNIFICANCE LEVEL 


SHORTEST SIGNIFICANT RANGE TABLE 
(IF A COMPARISON IS TO BE MADE ACROSS M (TRANSFORMED) MEANS, 
THE NECESSARY DIFFERENCE FOR SIGNIFICANCE IS S) 


M = 2 3 4 5 
S = 9.578 10.075 10.340 10.606 
YEAR RAW MEANS TRANS. MEANS NON-SIGNIFICANT GROUPINGS 

88 91.00 76.7647 
87 79.50 65 .2467 e 
86 68.50 56.2502 ee 
85 64.00 53.7166 oe 
84 49.00 44.4010 e 


DO YOU WANT TO MAKE ANOTHER RUN(YES OR NO)? 
=NO 
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EXHIBIT 17. 





SPSSX command file for two-way ANOVA with repeated measurements 


10$$S,J,T, MONI 

20$ : IDENT :SS01/FENNEC, $C210FKM 

30$ : SELECT : SPSS/X2.2/MAIN 

40$: PRMFL:11,R,S,SC325FM/RANGEMON/YYDATA 
SOFILE HANDLE DATAINPT /UNIT=11 

60DATA LIST FREE FILE=DATAINPT 

70 /TRANS FREQ1 TO FREQS 

75SET WIDTH=80 

80LIST VARIABLES=ALL 

85DO REPEAT FREQ=FREQ1 TO FREQS 
SOCOMPUTE FREQ=FREQ/20. 

1OOCOMPUTE FREQ=ARSIN (SORT (FREQ) ) 
LIOCOMPUTE FREQ=360.*FREQ/(2.*3.1416) 
11S5END REPEAT 

120MANOVA FREQ] TO FREQS/ 

122 WSFACTORS=YEAR(S) / 

124 WSDESIGN=YEAR/ 

126 ANALYSIS (REPEATED) / 

130 PRINT=<OMEANS (VARIABLES (FREQ1 TO FREQS) TABLES (CONSTANT) ) / 
140 PRINT=<SIGNIF (UNIV) / 

150 PRINT=PARAMETER (ESTIM) / © 
155 PRINT<TRANSFORM/ 

160 CONTRAST (YEAR) =POLYNOMIAL/ 

170 DESIGN/ 

180FINISH 

190$: ENDJOB 
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EXHIBIT 18. 





SPSSX results for two-way ANOVA with peated measurements 


e*eeereeree AR ALYSIS OF VARIANCE -- DESIGN zeeeeee 
COMBINED OBSERVED GRAND MEANS 
VARIABLE .. FREQ1 


WGT. 53.71648 
UNWGT . 53.71648 
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EXHIBIT 18 (CONT.) 


eee e*e2eCaRaALYsSsiIs OF VARIANCE -- DBSIGNR i***** * 
TESTS OF BETWEEN-SUBJECTS EFFECTS. 


TESTS OF SIGNIFICANCE FOR Tl USING UNIQUE SUMS OF SQUARES 


SOURCE OF VARIATION ss DF MS F SIG OF F 
WITHIN CELLS 1294.05 9 143.78 
CONSTANT 175680 .48 1 175680.48 1221.85 .060 


ESTIMATES FOR Tl 
--- INDIVIDUAL UNIVARIATE .9500 CONFIDENCE INTERVALS 


CONSTANT 
PARAMETER COEFF. STD. ERR. T-VALUE SIG. T LOWER -95% CL- UPPER 
1 132 .544515 3.79187 34.95491 -00000 123.96671 141.12232 











EXHIBIT 18 (CONT.) 





e*e*eees**ABALYStsS OP VARIABCS#S-- DESIGN iff ff #8 
TESTS INVOLVING ‘YEAR’ WITHIN-SUBJECT EFFECT. 


AVERAGED TESTS OF SIGNIFICANCE FOR FREQ USING UNIQUE SUMS OF SQUARES 


SOURCE OF VARIATION Ss DF MS F SIG OF F 
WITHIN CELLS 3954.31 36 109 .84 
YEAR 6028 .27 4 1507.07 13.72 .000 


ESTIMATES FOR T2 

--- INDIVIDUAL UNIVARIATE .9500 CONFIDENCE INTERVALS 

YEAR 

PARAMETER COEFF. STD. ERR. T- VALUE SIG. T LOWER -95% CL- UPPER 


1 -22.884498 3.70785 -6.17191 -00016 -31.27223 -14.49676 


ESTIMATES FOR T3 


PARAMETER COEFF. STD. ERR. T- VALUE SIG. T LOWER -95% CL- UPPER 


© 1 -5.6062508 2.58240 -2.17095 .05802 -11.44804 . 23554 


ESTIMATES FOR T4 


PARAMETER COEFF. STD. ERR. T- VALUE SIG. T LOWER -95% CL- UPPER 


i -1.7700235 2.67441 -.661864 52466 -7.81995 4.27990 


ESTIMATES FOR TS 
--- INDIVIDUAL UNIVARIATE .9500 CONFIDENCE INTERVALS 
YEAR 


PARAMETER COEFF. STD. ERR. T-VALUE SIG. T LOWER -95% CL- UPPER 


1 6.67559673 4.04567 1.65006 .13333 ~2.47633 15 .82753 














EXHIBIT 19. 





SPSSX command file for two-way ANOVA with more than one 
observation per block 


10$$S,J,T,MONI 

20$ : IDENT :SS01/FENNEC, SC210FKM 

30$ : SELECT :SPSS/X2.2/MAIN 

40S : PRMFL:11,R,S,SC325FM/RANGEMON/DDDATA 

SOFILE HANDLE DATAINPT /UNIT<11 

60DATA LIST FREE FILE=DATAINPT 

70 /TRANS YEAR REPL DENS 

74SET WIDTH<80 

7SLIST VARIABLES<ALL 

76COMPUTE DENS=<SQRT (DENS) 

120MANOVA DENS BY REPL(1,3) TRANS (1,2) YEAR(1,4)/ 
130 PRINT<OMEANS (TABLES (CONSTANT YEAR) ) / 

140 PRINT<SIGNIF (UNIV) / 

170 DESIGN«REPL WITHIN TRANS<1, TRANS VS 1, YEAR VS RESIDUAL, 
175 YEAR BY TRANS VS RESIDUAL/ 

180FINISH 

190$ : ENDJOB 











EXHIBIT 20. 


SPSSX results for two-way ANOVA with more than one observation per block 


eeee27eeQARaALYsSse Is opr VARIANCE -- DESIGN teeseee2s2 
COMBINED OBSERVED GRAND MEANS 
VARIABLE .. DENS 
HMEAN 
WGT. 4.45833 
UNWGT . 4.45833 


COMBINED OBSERVED MEANS FOR YEAR 


1 WGT. 3.83333 
UNWGT . 3.83333 
2 WGT . 2.50000 
UNWOT . 2.50000 
3 WG? . 6.16667 
UNWGT . 6.16667 
4 WGT. 5.33333 
UNWGT . 5.33333 


eeeceeeceaRALYseiItsSs opr VARIANC E -- DESIGN aeeeeee 
TESTS OF SIGNIFICANCE FOR DENS USING UNIQUE SUMS OF SQUARES 


SOURCE OF VARIATION ss DF MS F SIG OF F 

RESIDUAL 14.83 12 1.24 

YEAR 47.46 3 15.82 12.80 .000 

YEAR BY TRANS 7.46 3 2.49 2.01 .166 

ERROR 1 17.17 4 4.29 

TRANS 51.04 1 51.04 11.89 .026 
E-65 








EXHIBIT 21 


STPK results for two-way ANOVA with more than one observation per block 


STPK 


WHICH STAT GROUP WOULD YOU LIKE TO USE, 
TYPE A NUMBER FROM 1 TO 5, SOS, OR END. * 
=3 


TO CHANGE TO A DIFFERENT STAT GROUP, TYPE FINISH IN ANSWER TO 
-WHAT ANALYSIS DO YOU WISH TO PERFORM- 


TO TERMINATE STATPKG, TYPE END IN ANSWER TO 
-WHICH STAT GROUP WOULD YOU LIKE TO USE- 


IS YOUR DATA COMING FROM OTHER THAN TTY * 
=YES 

ENTER CATALOG FILE STRING FOR DATA FILE 
=SC325FM/RANGEMON/DDDATA 


SPECIFY THE NUMBER OF FACTORS * 
=3 


SPECIFY ONE-LETTER LABELS FOR EACH FACTOR IN THE FORM X/X/X/... 
=T/Y/R 


SPECIFY LEVELS (CATEGORIES) FOR EACH FACTOR IN THE FORM X/X/X/... 
=2/4/3 


SPECIFY INPUT CODE FOR YOUR DATA * 
=DF 


EACH INPUT RECORD SHOULD CONTAIN A SINGLE OBSERVATION PRECEDED BY ITS 
IDENTIFING SUBSCRIPTS - ONE FOR EACH FACTOR 


THE DATA WILL NOW BE READ FROM DISK IN A FREE FORMAT 


DO YOU WISH TO PRINT THE DATA JUST READ IN 
=YES 
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EXHIBIT 21 (CONT) 





FACTOR LEVELS OBSERVATION 
0. 

16 .00000 
0. 

4.00000 
25 .00000 
49 .00000 
9.00000 
64 .00000 
9.00000 
25 .00000 
1.00000 
16 .00000 
.00000 
36 .00000 
16.00000 
36 .00000 
16 .00000 
49 .00000 
9.00000 
25 .00000 
36 .00000 
64 .00000 
4.00000 
81.00000 


~ 
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wi 


DO YOU WISH TO MAKE ANY CORRECTIONS TO THE DATA * 
=NO 


DO YOU WISH TO PERFORM DATA TRANSFORMATION 
=YES 


SPECIFY TRANSFORMATION CODE * 
=SOS 


THE FOLLOWING CODES SIGNIFY TYPES OF TRANSFORMATIONS. . 
1 SQUARE ROOT OF X 


2 COMMON LOG OF X 

3 NATURAL LOG OF X 
4 EXPONENTIAL OF X 
9 x+cC (C IS A CONSTANT.) 
10 a?°C¢ 

11 1/X 


13 x ** Cc 
0 NO MORE TRANSFORMATION 
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EXHIBIT 21 (CONT.) © 


SPECIFY TRANSFORMATION CODE * 
=l 


SPECIFY TRANSFORMATION CODE * 


=0 
FACTOR LEVEL 
T 2 
Y 4 
5 3 
GRAND MEAN 4.45833 
LINE SOURCE OF SUM OF DEGREES OF MEAN 
NUMBER VARIATION SQUARES FREEDOM SQUARE 
1 T 51.04167 1 51.04167 
2 Y 47.45833 3 15 .81944 
3 TY 7.45833 3 2.48611 
4 Pp 14 .58333 2 7.29167 
5 TR 2.58333 2 1.29167 
6 YR 9.41667 6 1.56944 
7 TYR 5.41667 6 0.90278 
TOTAL 137.95833 23 


DO YOU WISH TO POOL COMPONENTS * 
=YES 


HOW MANY COMPONENTS IN POOL 1. TYPE O WHEN YOU WISH TO STOP * 
=2 


TYPE THE LINE NUMBERS OF THESE COMPONENTS IN THE FORM X/X/X/... 
=6/7 


LINE NO. 
8 POOL 1 14 .83333 12 1.23611 


HOW MANY COMPONENTS IN POOL 2 
=2 


TYPE THE LINE NUMBERS OF THESE COMPONENTS IN THE FORM X/X/X/... 
=4/5 


LINE NO. 
9 POOL 2 17.16€667 4 4.29167 


HOW MANY COMPONENTS IN POOL 3 








©} EXHIBIT 21 (CONT) 


DO YCU WISH TO COMPUTE F-VALUES * 
=YES 


TYPE LINE NUMBER OF MAIN EFFECT OR INTERACTION, AND LINE NUMBER OF 
ERROR TERM IN THE FORM XX/XX. TYPE 0/0 WHEN YOU WISH TO STOP 

=2/8 

F-VALUE..... 12.79775 (D.F.= 3, 12) 


TYPE LINE NUMBERS FOR NEXT F CALCULATION 
=3/8 


F-VALUE..... 2.01124 (D.F.= 3, 12) 


TYPE LINE NUMBERS FOR NEXT F CALCULATION 
=1/9 


P-VALUE..... 11.89320 (D.F.= i, 4) 


TYPE LIWE NUMBERS FOR NEXT F CALCULATION 
=0/0 


© DO YOU WISH TO POOL MORE COMPONENTS OR COMPUTE MORE F-VALUES 
=NO 


DO YOU HAVE OTHER ANALYSIS-OF-VARIANCE PROBLEMS TO SOLVE 


WHICH STAT GROUP WOULD YOU LIKE TO USE, 
OR END. * 
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EXHIBIT 22. © 


SPSSX command file for Chi-square test 


10$$S,J7,T,MONI 

20$ : IDENT :SS01/FENNEC, SC210FKM 

30$ : SELECT: SPSS/X2.2/MAIN 

40$ : PRMFL:11,R,S,SC325FM/RANGEMON /CHDATA 
SOFILE HANDLE DATAINPT /UNIT=11 

60DATA LIST FREE FILE=DATAINPT 

70 /YEAR CLASS COUNT 

7SSET WIDTH=80 

76LIST VARIABLES=ALL 

77WEIGHT BY COUNT 

80CROSSTABS VARIABLES=YEAR (84,86) CLASS(1,3) /TABLES=CLASS BY YEAR 
S9OSTATISTICS 1 

180FINISH 

190$ : ENDJOB 














EXHIBIT 23. 


SPSSX Results for Chi-square test 


---+-+++-+-++- CROSSTABULATION OF --------- 


CLASS 
BY YEAR 
cee ee 2.828 5 @ 2 6 es «8 2 6 2 2 24 6 2 2 & 2.86 8 PAGE 1OF 1 
YEAR 
COUNT I 
I ROW 
I TOTAL 
I 841 8sI 86I 
CLASS # # # <<<<=+-=-- etait  ieiieiieteatin selina . 
1 tI . 14 I io 6 31 
I I I I 29.0 
Hosen eee eo ++ eeee-  iieieiietete + 
2 I i4 iI 22 I is I 54 
I I I I 50.5 
Sitti Stites S leiietetetetie + 
3UCUd&g ' @ io 6I S$ 2 22 
I I I I 20.6 
Home e eee Hoo eeeeee Siete > 
COLUMN 28 46 33 107 
TOTAL 26.2 43.0 30.8 100.0 
CHI - SQUARE D.F SIGNIFICANCE MIN E.F CELLS WITH E.F.< 5 
1.14923 4 0.8864 5.757 NONE 
NUMBER OF MISSING OBSERVATIONS « 0 
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EXHIBIT 24. 





STPK results for Chi-square test 


WHICH STAT GROUP WOULD YOU LIKE TO USE, 
OR END. * 


TO CHANGE TO A DIFFERENT STAT GROUP, TYPE FINISH IN ANSWER TO 
-WHAT ANALYSIS DO YOU WISH TO PERFORM- 


-WHICH STAT GROUP WOULD YOU LIKE TO USE- 


TO TERMINATE STATPKG, TYPE END IN ANSWER TO 
IS YOUR DATA COMING FROM OTHER THAN TTY * 

| =YES 

ENTER CATALOG FILE STRING FOR DATA FILE 
=SC325FM/RANGEMON /ONECDATA 


DO YOU WISH TO LIST NAMES OF AVAILABLE ANALYSES * 


«YES 

THE FOLLOWING ANALYSES ARE AVAILABLE IN THIS GROUP.. © 
ep ADD, REPLACE, AND DELETE DATA 
TRANSFORMATION........... LOG, SQUARE ROOT, ETC. 

ELEMENTARY STATISTICS.... MEAN, STANDARD DEVIATION, ETC. 
CORRELATION 

CROSS TABULATION......... FREQUENCY TABLE FOR 2 CROSSED VARIABLES 
SCATTER DIAGRAM 

HISTOGRAM 

LINE PLOT................ VERTICAL LINE PLOT 

RANK CORRELATION 

CHI-SQUARE............... CONTINGENCY TABLE 

We ci ccubictecsceseds TWO SAMPLE MEANS 


WHAT ANALYSIS DO YOU WISH TO PERFORM * 
=CHI - SQUARE 


DO YOU WISH TO READ NEW DATA * 
«YES 


HOW MANY ROWS IN YOUR DATA MATRIX * 
=3 


HOW MANY COLUMNS * 
=3 











EXHIBIT 24 (CONT.) 





SPECIFY INPUT CODE FOR YOUR DATA * 
=DF 


THE DATA WILL NOW BE READ FROM DISK IN A FREE FORMAT 


DO YOU WISH TO PRINT THE DATA JUST READ IN 
=YES 
7.000 14.000 10.000 
14.000 22.000 18.000 
7.000 10.000 5.000 


DO YOU WISH TO CHANGE SOME VALUES * 
=NO 


DO YOU WISH TO PRINT EXPECTED FREQUENCIES 
=YES 


EXPECTED FREQUENCY IN EACH CELL 


WHAT ANALYSIS DO YOU WISH TO PERFORM * 
=FINISH 
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